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Manpower in Physics — A National Emergency 


N the second report! of the War Policy Com- 

mittee of the American Institute of Physics, 
issued July 8, 1942, emphasis is placed upon the 
alarmingly serious situation of national man- 
power in physics. The Institute’s Committee 
believes that this situation has already reached 
the proportions of a national emergency which 
may be summarized as follows: 


1. The design, production, operation, and maintenance 
of new physical instruments of warfare are essential to the 
successful prosecution of this war. (Examples: submarine 
and aircraft location, improvement of anti-aircraft fire, 
automatic fuses, signaling devices, magnetic mines, etc.) 

2. The need for more physicists for these purposes is 
large and urgent. 

3. The number of physicists in this country is small, 
only about 7000, and a substantial portion of these are 
already engaged in direct war work, leaving altogether too 
few to provide physics training for Army and Navy per- 
sonnel and those needed for war research and production. 

4. Training of physicists is not an easy or short-time 
process. 

5. Unless prompt, effective measures are taken the 
shortage of physicists will be disastrously acute and no 
adequate program for training new physicists can be 
effected. 


It is undoubtedly true that the Army, the 
Navy, and the War Manpower Commission also 
are beginning to recognize the existence of this 
emergency situation in physics. To meet it, the 
War Policy Committee is urging that these 
groups take the following steps: 





1 The first report of this committee was published in 
J. App. Phys. 13, 345-347 (1942). 


1. Arrange for teacher training to provide for the very 
great amount of physics teaching which will be needed, 
not only in producing physicists, but in connection with 
training programs of the Army and Navy in which it is 
essential to convey some knowledge of physics to over 
200,000 men and women within a year. 

2. Revise the situation of physicists with respect to 
Selective Service so as to assure students and teachers of 
physics the possibility of continuing their work without 
uncertainty. 

3. Provide loans, scholarships, or other assistance to 
well-qualified students who need aid to continue their 
training in physics. 

4. Any men in the Army and Navy with physics training 
who are not actually employing that training in their 
work should be transferred to positions where physicists 
are now urgently needed, either by detail or by discharge. 

5. Start a public relations program as to the meaning of 
physics and its importance in the war, this being necessary 
to secure public approval and understanding of the 
necessity of the preceding four steps. 


Already the teaching loads in universities have 
reached tremendous proportions. Example after 
example is being provided showing that present 
summer session enrollments in physics are three, 
four, and even five times the normal summer 
enrollment. Thus there is ample evidence of the 
need for action of the type which the Institute’s 
Committee has recommended. What the situa- 
tion in universities and colleges is going to be 
during the next academic year is difficult to 
predict. A statement of the teaching problems 
as they appear to the War Policy Committee is 
reprinted on page 477 of this issue. 
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Science In The U. S. S. R.* 


J. G. CROWTHER 





Mr. J. G. Crowther has an intimate knowledge of Soviet science and scientists, 
numbering many of the latter among his personal friends. He is thus particularly 
well qualified to write with authority upon the personalities, ideas, and achievements 
of the present-day successors of Lomonosov and Mendeléeff. The freshness and vigor 
of Russian thought are vividly illustrated by Rehbinder’s views that science should 
consciously aim at creating new industries and not be content with improving old ones. 





INCE the time of Peter the Great, there has 

always been notable scientific research in 
Russia. Peter founded the Academy of Sciences 
in 1724, after the French model. He and his 
successors depended a good deal for their 
importation from abroad. 
Daniel Bernoulli went to St. Petersburg in that 
year, and shortly afterwards he was followed by 
Leonhard Euler. The Academy from its founda- 
tion. had a very high standard, but from the 
nature of its origin it was an exotic intellectual 
ornament rather than a working organ of the 
state. In the last ten years especially it has been 
fundamentally reorganized. The principle of the 
Academy during the first two centuries of its 
existence was the appointment of men of ability, 
the more gifted the better, with handsome 
salaries and light duties. They were then 
expected to follow the promptings of their 
inspiration, and embellish history with new 
knowledge. 


academicians on 


The reorganized Academy has quite a different 
principle. It is intended to be a powerful organ 
in the advancement of human life. It has the 
duty of planning scientific activity in the 
U.S.S.R. so that it will achieve that end most 
effectively. The Academy was transferred to 
Moscow so that its directorate would have closer 
contact with the Government and its institutes 
would be in a safer strategical location. 

The idea of planning is to develop scientific 
research in a balanced manner, discover and 

* This article is being reprinted with the generous per- 
mission of Endeavour. The journal, Endeavour, is a new 
quarterly review, published in England by Imperial Chemi- 
cal Industries, whose stated purpose is to carry overseas 
news of the continuing vitality and progress of the sciences 


in these days when not only nations but the international- 
ism of the sciences are threatened. 
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encourage gifted individuals, and lay down lines 
of investigation consonant with the main 
interests of the population. For instance, plan- 
ning is intended to prevent overlapping and 
stampedes into fashionable subjects. It will see 
that there is division of attention 
between, say, nuclear physics and radiophysics. 
The periodical reports on the activities of all 
laboratories help to prevent good men from being 
overlooked. The coordination of the general lines 
of research with the main interests of the 
population taps a strong stream of intellectual 
energy that might have been dissipated in 
other directions. 

The Physico-Technical Institute directed by 
Professor A. F. Joffe at Leningrad is the most 
famous of its kind in the U.S.S.R. It is one of 
the group of institutes attached to the Academy 
of Sciences, and its activities are a part of the 
Academy’s general plan of research. It was 
started by Joffe in September, 1918 in a few 
rooms of the old Polytechnic Institute. It was 
moved in 1923 into a handsome building which 
had been started in 1914 as a home for retired 
engineers, but never completed. 

As the Institute grew, large blocks of new 
laboratories and machine shops were added to it. 
The original staff consisted of Joffe, Semenov, 
Tchernitchev, Dorfmann, Lukirsky, and Frenkel. 
The first aim of the Institute was the training of 
physicists who could become directors of new 
institutes and foci of new growths of science in 
distant parts of the country. As a result of this 
plan, physical institutes have been founded in 
Kharkov, Dnepropetrovsk, Tomsk, Sverdlovsk, 
and Samarkand. 

This development sprang from Joffe’s insti- 
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tute. He is accordingly admired as the founder 
of modern physics in the U.S.S.R. His sixtieth 
birthday was celebrated recently by a fine 
volume in the Academy of Sciences’ Journal of 
Physics, which consisted exclusively of a col- 
lection of papers by his former pupils. 

Joffe’s own researches in recent years have 
been on semiconductors. He regards these as 
substances in which electrons have a moderate 
degree of freedom, in contrast with insulators in 
which they have very little, and in metals in 
which they have a great deal. Their study 
throws light on the mechanism of the breakdown 
of insulators, because this happens too suddenly 
for direct investigation. On the other hand, it 
cannot be deduced from the properties of metals 
because the flow of current starts too easily. 
Another line pursued by Joffe is the strength of 
solids. It is found that the theoretical strength 
calculated from Born’s theory is several hundred 


times the observed strength. In his researches in 
this problem, Joffe and his colleagues observed 
that the tensile strength of a crystal of rocksalt 
was increased 20 times while immersed in hot 
water. He attributed this to the removal by 
solution of cracks on the surface of the crystal, 
which, according to the theory of Griffith, are a 
major cause of the observed weakness of glass 
and other solids. 

Stepanow’s suggestion that friction between 
adjacent planes in a crystal under strain leads to 
a momentary liquefaction has arisen out of 
Joffe’s work, and has helped to explain why a 
piece of metal may be suddenly twisted and 
remain strongly set in its new shape. The metal 
at the surfaces at the planes of slip is momen- 
tarily melted, and then sets quickly while the 
piece is in its new shape. 

The recent investigations of S. E. Bresler 
and P. A. Phinogenov on microscopic steel balls 
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Professor A. F. Joffe. 


for bearings is a characteristic product of the 
Joffe school. It follows from a formula for the 
deformation of spheres first given by Hertz that 
small balls may be as effective in bearings as 
large ones, because there are more of them. 
Accordingly, mass production methods of mak- 
ing microscopic steel balls were investigated. 
Soft carbon steel wire was blown through an 
electric arc into an atmosphere of nitrogen, 
producing a “‘fog’’ of steel particles in it. When 
examined, these were found to be shiny and 
from a few thousandths to a few tenths of a 
millimeter in diameter. The larger ones were 
prolate, and deviated by 5 to 10 percent from the 
spherical; while the smaller ones, owing to the 
more powerful effect of surface tension, were 
nearly perfect. Balls about 1/10 mm in diameter 
were separated by sieves. When mixed with oil 
they reduced friction 25 times as much as oil 
used alone. 

Perhaps the most remarkable achievement that 
has come from Joffe’s institute is Skobeltzyn’s 
discovery of the very swift particles in cosmic 
rays. In 1926 he began a study of the velocity 
of electrons ejected from radioactive sources by 
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photons, with the aid of a Wilson chamber 
surrounded by a magnetic field. He published 
photographs in 1927 of particles which shot 
through the field with such energy that their 
curvature could not be detected. 

The apparatus with which he secured this 
fundamental observation was of the most 
provisional character. The timing of the appa- 
ratus was controlled by an old Atwood’s machine, 
which seems to be about the only valuable 
application of this dreary invention on record. 

Skobeltzyn obtained several photographs of 
positrons, but was not lucky enough’ to be the 
first to recognize their nature. He has recently 
been working on the components of the cosmic 
radiation, and has argued against the theory of 
Heisenberg and Euler that the spontaneous 
disintegration of mesotrons is the mechanism 
responsible for the formation of the electronic 
component observed in the lower layers of the 
atmosphere, because, in his opinion, it fails to 
explain the observed divergence between the 
curves describing the increase in the intensities 
of the hard and the soft components with height. 

One of the most striking figures in Soviet 
science is N. N. Semenov. He is indeed one of 
the most remarkable men in the world, and is 
the chief founder of the contemporary theory of 
chain reactions. 

Semenov and Kapitza are old friends and 
fellow students. In the Leningrad flat of the late 
Mme. Kapitza, the scientist’s mother, there used 
to be a beautiful oil painting of the pair in 
discussion over an x-ray tube. It is to be hoped 
that this historic picture is safe. Semenov 
started research as a physicist, and studies of 
the mechanism of the breakdown of insulation, 
in which a progressive change proceeds through 
the material, led him to the study of chain 
reactions. The classical experiment which accel- 
erated his school on its brilliant career was made 
by Chariton and Walta in 1926. Berthollet had 
observed in 1797 that phosphorus and oxygen, 
in spite of their intense attraction for each 
other at ordinary pressures, do not combine if 
the pressure of the oxygen exceeds a certain 
higher limit. Joubert, in 1874, added that the 
reaction also stops if the pressure falls below a 
certain lower limit. Chariton and Walta found 
that if argon were introduced into the vessel, 
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the upper pressure limit was lowered. Semenov 
showed that the phenomenon could be explained 
by the theory of chain reactions. Many more 
discoveries in the same field have come from his 
institute. This was founded as a_ separate 
laboratory in 1931, an old secondary school 
building being adapted to the purpose. 

Chariton is an old Cavendish student, though 
his main achievements have been accomplished 
in chemical physics. He has recently published 
a characteristic paper on the inability of small 
charges of explosives to detonate. He made 
small calibrated glass tubes which were closed 
at one end, the other end opening out into a 
funnel about 1 cm in diameter at the top. 
These were filled with nitroglycerine, which was 
left in contact with a priming charge of lead 
azide at the open end. It was found that if the 
diameter of the tube was less than 2.08 mm the 
detonation did not pass from the wide funnel- 
shaped part of the tube into the narrower 
cylindrical part. With a 12-percent solution of 
methyl ether in nitroglycerine, the limit rose to 
about 3.6 mm. One may suppose that, below a 
certain diameter, the duration of the scattering 
of the outer part of the charge is less than the 
duration of the reaction, so that detonation does 
not occur. 

The Semenov school has contributed greatly 
to the strength of the Soviet Union, besides the 
advancement of knowledge, because its studies 
on chain reactions have provided experts on 
“knock” in engines, explosives, and other things 
of practical importance. 

Peter Kapitza is well known to many British 
scientists. The handsome Institute of Physical 
Problems has been specially built for him in 
Moscow. He has recently invented a turbine 
air-liquefier, in which the initial compression of 
the air is only 5 atmospheres, instead of 200 
atmospheres as in the Linde process. This is 
achieved by taking the energy out of the gas by 
making it do work in the turbine, instead of 
depending on the Joule-Thomson effect, which 
is of secondary magnitude. As the temperature 
of the air flowing through the turbine, unlike 
that of the stream in a steam turbine, is low, 
the working fluid is much more dense, being 
about five times that of steam at 250°C. Conse- 
quently, the turbine has to be designed to 
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Professor Peter Kapitza. 


withstand considerable centrifugal forces, like a 
water turbine. 

Air is compressed by a 50-hp compressor, and 
after passing through water and air coolers is fed 
at about 7 atmospheres into the turbine. This 
runs at 40,000 r.p.m., and develops about 4 hp. 
The pressure of the air drops 75 percent in its 
passage through the turbine, and if it enters at 
about —158°C it emerges at about —187°C, 
the boiling point of oxygen. 

The present apparatus produces about 30 kg 
of liquid air per hour, and with minor improve- 
ments should be as efficient as ordinary high 
pressure liquefiers. It is compact, owing to the 
absence of auxiliary equipment, and it delivers 
liquid air within the short period of twenty 
minutes from starting. 

Another interesting scientist is P. Rehbinder, 
who directs research on the chemistry and 
physics of surfaces in the Physical Institute of 
the Academy of Sciences in Moscow. He has a 
burning enthusiasm and energy for research, and 
original views on the role of science in ‘society. 
He considers that science should not be the 
handmaid of industry, but should consciously 
aim at the creation of new industries. These are 
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vastly more valuable than trifling improvements 
in conventional manufacturing processes. 

Rehbinder’s researches include striking in- 
vestigations on the phenomena of hardness. He 
has determined the hardness of specimens by 
deducing the strengths of the molecular fields at 
their surface from measurements of wettability, 
and he finds that his results agree with the usual 
methods of mechanical scratching. His method 
has the advantage of not damaging the specimen. 

He found that thin films are often much 
harder than solid blocks of the same material. 
For instance, a thin film of lead on glass may 
be twice as hard as glass itself, though block 
glass is thirty times as hard as block lead. In 
contrast with this, the application of other 
substances adsorbed by the surface may reduce 
the hardness. 

Rehbinder has recently been examining the 
effect of the addition of surface-active substances 
to the mechanical properties of soap curds. He 
has added homologues of a number of saturated 
alcohols, phenol, aniline, etc., to the curds of 
metallic soaps. He has found that the resistance 
of the curds to shear is reduced. He believes 
that the effect is partly analogous to the one he 
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Scientific education: a chemical laboratory at a university center in Daghestan. 








has already described, in which hardness of 
solids is reduced by adsorbed layers. He thinks 
that it should be possible to alter the mechanical 
properties of hard commercial soaps considerably 
by adding surface-active substances to them 
during manufacture and use. 

Perhaps the most gifted physicist in the 
U.S.S.R. is L. I. Mandelstamm, whose sixtieth 
birthday, like Joffe’s, was recently celebrated by 
a special volume in the Journal of Physics of 
papers by his pupils. He works in Moscow 
University and the Physical Institute of the 
Academy. Mandelstamm is modest, retiring, and 
self-critical. He does not like lecturing, and is 
cautious in publication. He discovered, in 
collaboration with Landsberg, the effect which 
now bears the name of Raman, in February, 
1928. While they were checking their results, 
Raman independently discovered the same effect 
and promptly published his observations in 
March, 1928. This is only part of Mandelstamm’s 
work. He has created a school of research in 
non-linear oscillations which is of the highest 
importance in radiophysics and aerodynamics, 
and helps to explain the strength of the Soviet 
Union in these fields. 
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The theory of non-linear oscillations was 
founded a century ago in connection with 
planetary perturbations, and was developed by 
Poincaré: Mandelstamm began the application 
of the theory of physical problems about fifteen 
years ago. Besides leading to improvements in 
radio circuits, it has inspired the invention of a 
new type of dynamo. It consists of a condenser 
which is placed in a circuit containing a self- 
inductance. The plates of the condenser are 
shaken by a mechanical motor, so that the 


distance between them is always changing. 
This varies the capacity of the condenser and 
produces a current in the circuit whose period is 
one-half that of the mechanical oscillator. Thus 
the energy of the mechanical oscillator is 
converted into a high voltage electric current. 

Such are a few of the Soviet scientists and 
their researches. There are hundreds of others, 
in very different branches of science. They are 
an example of the quality and energy of the 
scientific workers in the U.S.S.R. 








An Estimate of the Teaching Problems in Colleges and Universities for 
1942-43 as Made by the War Policy Committee of 
the American Institute of Physics 


(Report of July 8, 1942) 


66 HE Committee has studied plans of the Army and the Navy to enlist a large 

fraction of next year’s college and university students in training programs 
to be carried on at the institutions. The programs include the Army Enlisted 
Reserve and the Navy V-1 program. The Services rely on these programs to provide 
large pools of officer material and can be expected to make every effort to recruit 
them to the necessarily large enrollments to meet the essential needs. The Navy 
program requires a substantial course in physics; the Army program is such as to 
increase above normal the number of students electing physics. On the basis of 


joint announcements of the Services and conferences with their representatives, the 


Committee has estimated that the teaching load in college grade physics will in 
1942-43 be from two to three times the highest ever before sustained. This teaching 
load will fall on faculties necessarily depleted for war research and, in some in- 
stances, wastefully dissipated through the direct action or indirect influence of the 
Selective Service Act and the lack of knowledge in some local boards as to the im- 
portance of physics and its wartime role. 

“The Committee welcomes the veryspecific directives which have been issued by 
Selective Service Headquarters, but urges that high officials of the Army and the 
Navy issue statements emphasizing the necessity for students and teachers to 
continue their present work and affirm that in so doing the individuals concerned 
are performing the highest and most patriotic duty open to them. Such public 
statements, as well as occupational deferment policies, should particularly include 
graduate students, since these students are already contributing heavily to the 
teaching of physics and are being drawn more and more into war research. 

‘Whatever actions may be taken by Federal authorities (and some are known to 
be in preparation) designed to lessen the prospective shortage of teachers, the 
Committee is convinced that they cannot be wholly adequate. The Committee, 
therefore, urges college and university administrations and faculties to use every 
expedient they can devise to prepare for the coming flood of enrollments in physics 
courses. This advice is all the more urgent in view of the fact that the war research 
programs must continue to expand and that additional faculty members will have 
to be called away from their teaching duties.”’ 
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Response by Dr. William D. Coolidge in Accepting the Duddell Medal 
at the Meeting of the American Physical Society, May 1, 1942 


YOUR EXCELLENCY, MR. PRESIDENT, MEMBERS 
OF THE AMERICAN PHYSICAL SOCIETY, LADIES 
AND GENTLEMEN: 


AY I first thank Your Excellency for your 

great courtesy in coming here and giving 
me the pleasure of receiving this honor at your 
distinguished hands. This is generous compensa- 
tion for the disappointment which I naturally 
felt because present conditions made it inadvis- 
able for me to attend a meeting of The Physical 
Society in London to receive the award directly 
from their president. 

I am deeply conscious and appreciative of the 
honor of being made a Duddell Medalist, and for 
several reasons. The medal takes high luster from 
the great scientific achievements of Duddell, in 
whose honor The Physical Society established the 
memorial fund—a luster enhanced by the emi- 
nence of that society itself and by the conspicu- 
ous merit of former recipients. It carries special 
significance and value for me in coming from 
England, the home of my ancestors, the source 
of many great achievements in science, and the 
stalwart guardian of what are to me the finer 
things of life. 

To return for a moment to the work of Dud- 
dell. Had he done nothing more than develop the 
magnetic oscillograph, his contribution to the 
electrical science and art would have been out- 
standing. One cannot visit a laboratory today 
where electrical researches or electrical engineer- 
ing developments are in progress without being 
greatly impressed by the importance of the role 
played by the oscillograph which yields so much 
essential information obtainable in no other way. 
For frequencies within its range, the magnetic 
oscillograph, substantially as developed by Dud- 
dell, is still in general use. For higher frequencies, 
the cathode-ray oscillograph is used, but here 
also credit must be given to Duddell, for his mag- 
netic oscillograph served to blaze the trail. The 
art of radio also is heavily in’debt to Duddell for 
his discovery of the singing arc, his design of a 
high frequency alternator, and of his thermo- 
ammeter, all of which contributed to the early 
development of wireless telegraphy. 
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To make clear the limitations of my part in the 
work on which the award was based, I must 
briefly outline the conditions under which the 
work was done. 

Our laboratory, at the time the ductile tung- 
sten work was started, was only five years old. 
This laboratory sprang from the vision of Dr. E. 
W. Rice and his associates. He was then in charge 
of General Electric engineering and had seen 
clearly the extent to which the products of our 
factories had been derived from the results of 
scientific research, such as had been carried out 
notably by English physicists. 

The laboratory was started in the hope that it 
might contribute new facts and principles on 
which new products could be based. It was one 
of the earliest industrial research laboratories in 
this country and was established primarily for 
fundamental research. It was itself an experi- 
ment, but an experiment carried out under ex- 
tremely favorable conditions, as both its sponsors 
and its director, Dr. W. R. Whitney, were men 
of great courage, optimism, and faith in the 
power of scientific research methods. The strong- 
est single psychological factor conducing to the 
success of the laboratory was the spirit of co- 
operation which Dr. Whitney was able to build 
up not only within the laboratory but also 
between the staff and the engineers of the 
Company. 

That cooperation was so close and so generous 
and so effective that a very large share of the 
credit for the ductile tungsten and x-ray work 
belongs to those who were associated with me in 
that work. Credit is also due to the courageous 
and liberal management of the Company for giv- 
ing its continued support to such work as the 
ductile tungsten effort, which over a long period 
of years was exceedingly expensive and of ques- 
tionable promise. 

While the factory process for making ductile 
tungsten is today fundamentally the same as that 
developed in the laboratory, it has undergone 
constant improvement and been given increased 
utility by our engineers for the various applica- 
tions which have been found for the product. 
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Although this process did not constitute the 
first application of powder metallurgy, it seems 
to have been at the time the most important one 
which had been made, and, for this reason it'may 
perhaps be fairly credited with having been in no 
small measure responsible for the subsequent de- 
velopment of the art of powder metallurgy now 
widely applied to other metals, and to mixtures 
of metals and refractory compounds. 

My participation in the work on x-rays, which 
forms the other basis for the award, followed as a 
result of the keen interest which I had felt in 
X-rays ever since the announcement of their dis- 
covery, also from the work which we did in 
adapting wrought tungsten to use as a target in 
the earlier x-ray tube in place of platinum, and 
from Dr. Irving Langmuir’s studies on the 
thermionic emission of electrons from tungsten 
filaments in high vacuum. I had the advantage 
also of intimate firsthand knowledge of the limi- 
tations of the x-ray tube of that time and of the 
source of those limitations. Stable and control- 
lable thermionic electron emission rendered un- 
necessary the gas content which alone had been 
responsible for the most troublesome limitations 
of the earlier tube. At the same time it provided 
the much desired flexibility of control and opera- 
tion. I was also fortunate in having readily avail- 
able not only tungsten wire for the cathode 
filament but also large masses of wrought tung- 
sten which were at the time indispensable for the 
target. For when our first successful tubes were 
made, we did not know that it would ever be 
possible to degas sufficiently a composite target 
of copper and tungsten such as had been generally 
used in the earlier type of tube. 

I cannot forego this opportunity to pay my 
tribute to Sir Owen Richardson who, more than 
any other, had laid the groundwork for our 
knowledge of thermionic emission. 

From the beginning, it was clear that the high 
voltage generating equipment which had been 
developed for use with the earlier tube was ill- 
suited to the hot cathode type and for this reason 
later developments of the hot cathode tube and 
of the high voltage source have, even up to the 
present time, gone hand in hand. 

While some hot cathode tubes of the first so- 
called universal type are still in use, this design 
has generally been supplanted by special tubes 
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giving greatly improved performance for various 
specific uses. 

For example, types have been developed for 
operation on unrectified alternating current. The 
characteristics of the tube have permitted vari- 
ous forms of equipment in which the tube and the 
entire high voltage circuit are completely en- 
closed within grounded metal, thus eliminating 
all danger of electric shock and facilitating the 
attainment of adequate x-ray protection. For 
radiography, definition has been greatly im- 
proved in various ways, notably by the use of 
elongated focal spots and more recently, by this 
combined with rapid rotation of the target. 

For medical therapy and for industrial radi- 
ography, tubes for higher voltages than could be 
used with the first design have been developed, 
first by merely increasing dimensions. Then, be- 
cause of troublesome field currents arising from 
too high potential gradients at the cathode, we 
were led to depart from the simple two-electrode 
structure and to use a multiplicity of tubular 
accelerating electrodes in addition to cathode and 
anode. This apparently makes it possible to pro- 
vide tubes for voltages as high as can be produced. 

A recent development in the generation of high 
voltage to supply such tubes is the resonance 
transformer. This has no iron within the primary 
and secondary coils, thus permitting the x-ray 
tube to be located here, where it is electrostati- 
cally shielded by the surrounding coils. The size 
and weight of such equipment have been greatly 
reduced by substituting for oil, a compressed gas 
having high dielectric strength. 

For x-ray diffraction work, the recent develop- 
ment of tubes with linear focus combined with 
metallic beryllium windows results on the one 
hand in great reduction of the required exposure 
time, and on the other in a considerable exten- 
sion, in the long wave-length direction of the 
available spectrum. 

It is exceedingly pleasant to think of any useful 
purpose subserved by war. World War No. I did 
greatly accelerate both the development and the 
introduction to use of the hot cathode tube. The 
present war has already led to the extensive use 
of the photoroentgenographic method for the 
chest examination of recruits for the armed forces 
and this, because of its greatly reduced cost, ex- 
cellent results, and convenience, is certain to find 
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widespread use after the emergency. The present 
war is also stimulating enormously the applica- 
tion of x-rays to industrial uses, both radio- 
graphic and fluoroscopic. As an illustration of 
this, although the first portable million-volt 
radiographic x-ray outfit went into use only about 
a year ago, the fortieth one will soon have been 
completed. These million-volt outfits make pos- 
sible the rapid examination of even as difficult 
subjects as steel castings and welds having a 
thickness of as much as eight inches, and thereby 
are contributing most helpfully to many indus- 
trial war activities. As another example, fluoro- 
scopic outfits are being built and used for the 
rapid examination of small aluminum die cast- 
ings, thus permitting immediate rejection of 
many whose defects could otherwise be seen only 
after a considerable amount of machine work had 
been performed on them. 

To me, Your Excellency, the greatest value of 
this award lies in the fact that it is one of the 
many evidences of the cordial and friendly spirit 
which has for years existed between British and 
American scientists. This has been evidenced 
anew by the close and effective cooperation in the 
present war work between the scientists of the 
British Commonwealth and those of this country. 
It has been greatly facilitated through your vari- 
ous war agencies in England, the National Re- 
search Council of Canada, and our Office of Scien- 
tific Research and Development. Your scientists 
are bringing us full reports of the vast amount of 
splendid work which you are doing in support of 
the war effort. We are making good use of it and 





are in turn keeping you informed concerning our 
work. This augurs well, both for the emergency 
and for the more leisurely pursuits to follow. 
This mutual aid in our common Cause is strength- 
ening the bonds between our nations, and we may 
hope that, when victory is won, those bonds may 
serve to hold us together in the task of rehabili- 
tating the world on a firm basis of liberty, 
security, and opportunity for all. 

In every country of the world there are men 
of good will, who desire peace and friendly rela- 
tions with other nations. Misrepresentative gov- 
ernments may gain temporary ascendancy, and 
by ruthless aggression and acts of cruelty bring 
upon themselves the hatred of their victims, but 
there remains, though temporarily submerged 
and silenced, the peace-loving mass of human 
beings who comprise the bulk of the nation, and 
who will respond in kind to an approach in the 
spirit of friendship and mutual helpfulness. 

It is a far cry today, but we must certainly look 
hopefully to the time, and labor diligently for it, 
when such relations of mutual esteem and friend- 
ship as now exist not only between the scientists 
but also between the peoples of the British 
Commonwealth of Nations and the United States 
of America can spread to all peoples and to all 
countries. 

In closing, Mr. President, I wish to express to 
the American Physical Society my gratitude for 
graciously providing this very pleasant setting 
for this occasion and thus lending indorsement 
to the honor conferred upon me by The Physical 
Society of London. 


ESEA RCH, after all, is nothing but prospecting for knowledge. Whenever 
someone sets out on an organized program to uncover new knowledge about 





a subject of which he knows very little, he is engaging in research. There is nothing 
mysterious about it, nothing highbrow. It is the most natural thing in the world. 


CHARLES F. KETTERING in “Prospecting for Knowledge” 
in the Bakelite Review, January, 1942 
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Résumé of Recent Research 











A New Type of The optical system of 
Variable Focus Lens* the eye is made up of 

two main elements: a 
fixed focus lens (the cornea) and an auxiliary lens 
of variable focus (the so-called ‘‘crystalline”’ 
lens). While great progress has been made in the 
design of fixed-focus lenses to aid defective vision, 
a practical lens of variable focus has not been 
made available. Attempts to produce lens varia- 
bility have resulted either in complex optical 
systems, such as the ‘“‘zooming focus’’ devices 
and the new oscillating lens systems used on 
cinema cameras, or in variable but optically 
unsatisfactory single lenses. The single lenses, for 
the most part, consist of liquid bounded by two 
transparent walls which are segments of spheres 
placed back to back. By constricting the pe- 
riphery of these liquid-filled semi-spheres, or by 
increasing the pressure of the liquid within them, 
the attempt is made to maintain spherical sur- 
faces while changing the radii thereof. These two 
things cannot simultaneously be done without 
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Fic. 1. 


* Patent No. 2,269,905 (1942). 


VOLUME 13, AUGUST, 1942 





actually altering the surface areas. In the absence 
of suitably elastic and transparent material such 
attempts result in distortion of the lens surfaces. 











Fic. 2. 


It is possible, however, to change the radius 
of a cylindrical segment without appreciably 
affecting its area and to produce a spherical 
refracting effect by means of two cylindrical 
segments, crossed at right angles. A lens embody- 
ing these two principles was designed and con- 
structed at Ohio State University.! 

The essential refracting unit consists of two 
plane disks of glass with clear liquid held between 
them by a rubbe: dam which connects the edges 
of the glass surfaces. For flexibility the optimum 
thickness of the glass was found to be about 0.07 
mm. Clamped only by rubber and protected by 
glass cover plates, the two thin disks are not 
easily broken. 

To select any dioptral power within the limits 
of the mechanism (+3.50 diopters in the first 
instrument), the circular plate of Fig. 1 is 
rotated until the desired figure is opposite a 
fixed indicator arrow. The reverse side of the 
device is identical to the side shown except that 
this second side is oriented at right angles to 
the first, so that the axes of the two variable 
cylinders will be crossed. By suitable rotation of 
the two circular plates with respect to each 
other, cylindrical, spherical, and sphero-cylin- 
drical powers may be obtained. 

This instrument (Fig. 2) was designed pri- 
marily for use in measuring errors of refraction 
of the eye to replace the ordinary method of 


1 Robert Graham, J. Opt. Soc. Am. 30, 560-563 (1940). 
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measuring such errors by trying lens after lens 
before the eye until a combination is found which 
neutralizes them. In orthoptic or eye-training 
devices rotary prisms are frequently employed to 
alter gradually the amount of convergence of the 
two eyes while the accommodation remains con- 
stant. The variable lens now complements this 
practice by making it possible to induce gradual 
changes in the accommodation while the con- 
vergence remains constant. The lens may also 
find uses in certain types of cameras, where the 
proper focus may be procured merely by altering 
the focal length of the lens instead of by changing 
the separation between lens and film, as is done 
at present. 

It also seems possible that variable lenses as 
described above could be made small enough so 
that, in certain cases, they could be worn before 
the eyes to replace bifocal lenses. 


New Books 








Photoelasticity, Volume I 


By Max Mark Frocut. Pp. 411+xxvi, 23}16 cm. 
John Wiley and Sons, Inc., New York, 1941. Price 
$6.00. 

This is the first of two volumes on the subject of stress- 
optics. It covers the “essential, basic and _practical’’ 
aspects of the subject and is intended for engineering 
reference and instruction. More advanced topics have been 
reserved for Volume II which is in preparation. 

The photoelastic technique has progressed materially 
since Coker and Filon published their treatise. No other 
text has been available in English, although there is a 
voluminous but somewhat scattered literature on the 
subject. The present work therefore fills a need which has 
often been expressed but which has not been filled because 
of the labor involved in organizing what is essentially the 
first text to place stress-optics on an engineering basis. 

The treatment includes a large number of illustrated 
examples of photoelastic analysis. This tends to clarify 
material which deals with basic principles not currently 
well understood by the engineering profession. The author 
points out that conventional analysis often fails to provide 
detailed stress distributions in machine and _ structural 
parts and that in many such cases the photoelastic pro- 
cedure offers the only reasonable approach. Many photo- 
graphs of stress patterns are included. These exhibit a 
technique and perfection for which the author is well 
known among his associates. 

A significant lack in Volume I is an adequate treatment 
of the conditions of elastic similitude and it is hoped that 
this will be presented in some detail in Volume II. In other 
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respects however, the reviewer wishes to indicate great 
admiration for the excellent manner in which the subject 
has been treated. 


R. WELLER 


Naval Ordnance Laboratory 


Acoustic Design Charts 


By FRANK Massa. Pp. 228+xiv, Charts 107, 23316 
cm. The Blakiston Company, Philadelphia, 1942. 
Price $4.00. 

This book, in a convenient size for reference purposes, 
contains 107 charts of useful data in the field of acoustical 
engineering. These are divided into ten sections according 
to subject matter. Each chart is explained by a short 
paragraph of text accompanying it, and its application is 
illustrated by a sample problem with solution. The book 
cannot claim to present any original material that cannot 
be found elsewhere in standard textbooks and in the 
literature. However, these charts are extremely handy for 
the person working in the acoustical field, who probably 
has drawn up similar charts in connection with his work 
—but not as many, nor as well designed and conveniently 
arranged as in this book. It is possible to think of addi- 
tional charts that might be included, for instance direc- 
tional radiation characteristics are given for a vibrating 
piston and ‘a thin ring, but the well-known line array is 
omitted. The comment could also be made that no deriva- 
tions of the charts are given and that charts based on 
experimental data are included without distinction from 
those representing exact mathematical relations. 

These are very minor comments which do not detract 
from the value of the book as a whole, which should be a 
welcome addition to the library of all those interested in 
the design and construction of electro-acoustic apparatus. 


E. DIETZE 
Division of National Defense 


Columbia University 


Symposium on New Methods for Particle Size 
Determination in the Subsieve Range 


Pp. 111+Vv, Figs. 25, 15X23 cm. Published by the 
American Society for Testing Materials, Philadelphia, 
1941. Price $1.25. 

The papers and discussion presented at the Symposium 
on New Methods for Particle Size Determination in the 
Subsieve Range, held during the 1941 Spring Meeting of 
the American Society for Testing Materials, have recently 
been made available in book form. The titles of the papers 
and their authors are as follows: Methods for determining 
particle size distribution by Herbert E. Schweyer and 
Lincoln T. Work; Shape and surface of fine powders by 
the permeability method by P. C. Carman; Surface area 
of Portland cement by Paul S. Roller and P. V. Roundy, 
Jr.; An improved hydrometer method for use in fineness 
determinations by Alexander Klein; The determination of 
subsieve particle size distributions by sedimentation 
methods by S. W. Martin; The electron microscope in the 
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determination of particle size characteristics by James 


Hillier; A new method for measuring the surface areas of 
finely divided materials and for determining the size of 
particles by P. H. Emmett; Specific surface of pigments 
by adsorption from solution by Warren W. Ewing. 

Copies may be obtained from A.S.T.M. Headquarters, 
260 S. Broad Street, Philadelphia, Pennsylvania. 


Here and There 


Honors and Honored 








Professor Albert Ball, Chairman of the Department of 
Physics at Cooper Union, who is retiring after 37 years of 
service, Was given a testimonial dinner on May 28 by past 
and present members of the faculty of that school. 

Professor Frank M. Simpson, Chairman of the Depart- 
ment of Physics at Bucknell University, retired on July 1 
with the title Emeritus after 35 years of service in these 
capacities. ’ 

Dr. Richard C. Tolman, Professor of Physical Chemistry 
and Mathematical Physics at the California Institute of 
Technology, and Vice Chairman of the National Defense 
Research Committee, received the honorary degree of 
Doctor of Science at the commencement exercises of 
Prnceton University. 


* 
What Ho, the Attic! 


Research workers seeking instruments-required in their 
work but difficult to find are invited to communicate with 
D. H. Killeffer, Chairman of the newly appointed Com- 
mittee on the Location of New and Rare Instruments of 
the National Research Council. The plan of the Com- 
mittee’s activity is to assist in locating needed instruments 
of types not ordinarily available through accustomed 
channels. 

Assistance is particularly desired from owners and 
builders of instruments falling within the new or rare 
categories which might be made available to others through 
sale or for temporary use under mutually satisfactory 
conditions. 

Communications should be. addressed to D. H. KIL- 
LEFFER, 60 East 42nd Street, New York, New York. 


* 
New A.S.T.M. Officers 


Recently elected officers of the American Society for 
Testing Materials for 1942-1943 are as follows: President 
(one-year term), Professor H. J. Ball, Lowell Textile 
Institute; Vice President (two-year term), Mr. P. H. 
Bates, National Bureau of Standards; Members of the 
Executive Committee (two-year term), Dr. R. P. Anderson, 
American Petroleum Institute; Dr. Maurice H. Bigelow, 
Plaskon Company, Inc.; Mr. J. H. Foote, The Common- 
wealth and Southern Corporation; Mr. Alexander Foster, 
Jr., Warner Company; Mr. Lawford H. Fry, Edgewater 
Steel Company. 
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Fellowships Available 


Two fellowships are available in the Department of 
Geology, Massachusetts Institute of Technology, in the 
fields of Mineralogy and Structural Crystallography. One 
of these, a teaching fellowship, carries a stipend of $500 
per academic year plus free tuition. Up to eighteen hours 
per week may be required of the holder of this fellowship. 
The duties consist of teaching and other assistance in 
connection with laboratory instruction in mineralogy and 
possibly optical crystallography. Preference will be given 
to a student specializing in mineralogy or crystallography. 

Another fellowship, concerned chiefly with carrying out 
investigations in structural crystallography and _ allied 
fields, requires an investigator with a fairly good back- 
ground in crystallography and preferably having had some 
experience in the x-ray study of crystal structures. This 
fellowship, which is of a more permanent nature, requires 
the full time of the incumbent and carries a stipend of 
about $2000. 

Anyone interested in either of these fellowships should 
communicate with Professor M. J. Buerger, Department 
of Geology, Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts, giving full details regarding their 
qualifications. 


* 


Institute of Physics, London 


During the past year the Institute of Physics in London 
has established an Industrial Radiology Group and an 
Electronics Group, the first such subject groups to be 
formed under the auspices of the Institute. 

Recently elected officers who will assume duties on 
October 1, 1942 include: President, Sir Lawrence Bragg; 
Vice Presidents, Dr. W. Makower and Mr. T. Smith; 
Honorary Treasurer, Major C. E. S. Phillips; Honorary 
Secretary, Professor J. A. Crowther; Ordinary Members of 
the Board, Professor J. Chadwick, Professor J. D. Cock- 
croft, Mr. D. C. Gall, and Mr. E. B. Wedmore. 


* 


Graduate Courses in Ultra-High Frequencies 


The Graduate Electrical Engineering Department of the 
Polytechnic Institute of Brooklyn has again offered a 
group of five courses under the E.S.M.D.T. program for 
graduate students of physics or electrical engineering, as 
follows: Measurements at Ultra-High Frequencies; Cath- 
ode-Ray Oscillography; Generation and Reception of 
Ultra-High Frequencies; Ultra-High Frequency Antenna 
Systems and Reflectors; Ultra-High Frequency Wave 
Guides and Horns. 

The program, which started the week of June 15, extends 
over eight weeks, each course meeting once a week for 
eight three-hour sessions and once a week for laboratory 
work of from five to eight sessions. 
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Theory of Vehicle Suspensions 


PIERRE ERNEST MERCIER 
Electro- Mechanical Research, Inc., Houston, Texas 
(Received April 20, 1942) 


The present theory of the suspension of vehicles is established on the assumption of perfect 


symmetry and of a linear dependence between elastic forces and the displacements of the 


wheels and deflections of the tires. The path of the vehicle rolling on uneven ground is supposed 
nearly straight, but no other important limitations are made in the fundamental hypothesis, 
so that the results cover any kind of suspension. The relative magnitude of the seven vibrations 
corresponding to bobbing, rolling, and pitching motions of the body, and to the four specific 
vibrations of the wheels, are discussed in view of their combined effects on the displacements 


of the body of the car. The general conclusion, which attempts to separate as much as possible 


the respective frequencies to avoid beats and also directly to improve comfort, is in accord 
with several practical experiments made in the industrial field. 


INTRODUCTION 


HE object of this paper is to present a 

theory of the suspension of vehicles, i.e., 
to analyze the motion of the vehicle body and 
wheels caused by translation over uneven ground. 
Conclusions are drawn from this analysis which 
are of considerable interest in practical suspen- 
sion design. All such conclusions are in agreement 
with experimental evidence. 

Passage over uneven ground sets up oscilla- 
tions both of the wheels and of the sprung mass. 
The former occur at much higher frequencies, 
but nevertheless have an important effect on the 
latter; this may seem paradoxical at first, but 
becomes clear as the theory is developed, and is 
in agreement with experimental data. 

The sprung mass in general has three degrees 
of freedom, corresponding to motions in vertical 
translation or “bounce,” pitch, and roll, respec- 
tively. Particular emphasis is laid on the im- 
portance of the these 
modes. 


relationship between 


ASSUMPTIONS AND BASIC DATA 
Springing and Damping 


A sprung vehicle consists of a body supported 
on wheels to which it is connected by elastic 
suspension members, the wheels being fitted 
with tires which are also capable of elastic 
deflection. 
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The elastic suspension members (which may 
be leaf or coil springs, torsion bars, or air 
cylinders) are generally supplemented by some 
form of shock absorber or ‘‘damper.’’ These are 
not usually made so powerful as to suppress 
oscillations altogether, as this would have a 
detrimental effect on comfort. 

Damping forces of reasonable magnitude have 
but little effect on the frequency of vibrations, 
and in this particular problem even the effect on 
amplitude is comparatively small. For the sake 
of mathematical simplicity damping forces will 
therefore be left out of account. This simplifica- 
tion is further justified by the fact that when 
elastic characteristics are rationally determined 
the required degree of damping is quite low. 
The general validity of the conclusions will not 
be affected by this simplification. All loads are 
assumed directly proportional to displacements. 

The effective mass of each wheel is equal to 
the mass of the wheel and a fictitious mass which 
represents the inertia of the wheel carrying 
members; secondary reactions due to unbalance 
of the members are not taken into consideration. 
This assumption is quite evident for independent 
wheel springing suspension and admissible for 
orthodox suspension. 

But we shall deal with the most general type 
of suspension of which the orthodox axle sus- 
pension and the independent wheel springing 
are but particular cases, as will appear presently. 
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Type of Suspension 


The type of suspension is determined by the 
interaction between different wheels; i.e., the 
load in the suspension member corresponding to 
one wheel caused by the position and displace- 
ment of other wheels. In an independent sus- 
pension, such interaction is entirely absent if we 
suppose the spring mass non-affected by the 
motion of the wheels: that means that in an 
independent suspension there is no direct inter- 
action. In the axle type of suspension there is 
an interaction between the two front wheels and 
also between the two rear wheels. In this case, 
the principal parts of such interaction are elastic 
interactions due to the ratio between the mutual 
distance of both springs (corresponding to one 
axle) and the track itself (mutual distance of the 
two wheels on the same axle) and the secondary 
reactions mentioned above are not’ very 
important. 

One might, however, also conceive suspension 
giving interaction between the two wheels on 
the same side, or any other combination. The 
most general suspension will have interaction 
between all wheels; such suspensions can be 
realized in practice and more will be said about 
them later. For example, mechanical means can 
be developed in such case by connecting members 
working preferably in pure tension or by using 
torsional members. To reduce the weight and 
the dimensions of such systems the magnitude 
of the displacements communicated to the 
connecting members by the wheel motions are 
generally small in comparison to the wheel 
displacements themselves; so the specific differ- 
ences between the several types of inertia 
effects involved in rolling, pitching, etc., are 
very small despite the stress applied to those 
connective members, and the assumption that 
secondary reactions are negligible is also true in 
this case if the fictitious mass added to the mass 
of each wheel is correctly adjusted. 


Four-Wheel Vehicles 


The treatment will be restricted to four-wheel 
vehicles. For the sake of simplicity, certain 
symmetry conditions are assumed, viz: 

(a) Tracks of front and rear wheels are equal. 

(b) All wheels and tires are identical. 
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(c) The vehicle has two vertical planes of 
symmetry (fore-and-aft and transverse). 


The last assumption means that the c.g. is 
taken to be midway between front and rear 
wheels, and that front and rear suspension 
members are the same. Although this may seem 
an appreciable departure from reality (as far as 
practical vehicles are concerned), any more 
general assumption would lead to complications 
which would obscure the results without, how- 
ever, modifying their meaning to any appreciable 
extent. Furthermore, the general method given 
here may be extended if desired to cover any 
special case. 


DEVELOPMENT OF THEORY 
Basic Displacements 


Wheels are numbered in rotation round the 
vehicle: 1, 2, 3, 4, and corresponding suffixes 
will refer to the appropriate wheel, suffix ‘‘0”’ 
referring to the complete vehicle. 

Figure 1 shows the relative positions and 


x 





Zz 


Fic. 1. Axes of reference. 


displacements of the wheels and the axes of 
reference. Rotations corresponding to roll and 
pitch are denoted by ® and @, respectively. It 
will be noted that the positive direction of @ is 
the opposite of that usually taken (in order to 
maintain symmetry in the expressions for roll 
and pitch). 

The plane of reference (XY plane) is the 
ground when even. Unevenness of the ground 
may be denoted by rises or “motions” under 
the wheels, the motions under the four wheels 
being denoted by 21, 22, 23, and 24, respectively, 
the mean rise being denoted by 2p. 
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Let us consider any one wheel, say wheel 1. 
(See Fig. 2.) When the ground under this wheel 
rises by an amount —2, this is taken up by 
compression of the tire, deformation of the 
suspension member, and rise of the sprung mass 
above the plane of reference. Let these be 
denoted by —2;', —2)"’, and —2z,’", respectively, 
for the wheel 1 (similarly, with suffixes 2, 3, and 4 
for the other wheels). We must then have 


ay =a +21’ +21”, (1) 


and similarly for the other wheels (the quantity 
z;"’ is the deformation of the suspension member 
corresponding to wheel 1 as measured at the 
wheel center, i.e., the vertical movement of the 
wheel center relative to the sprung mass). 

All quantities are measured from the static 
position as origin. Similarly, loads are measured 
from the static loads as origin, since static 
deflections and gravity loads have no effect on 
vibration phenomena of the type considered, 
which are only concerned with displacements 
from the static position. It should be noted that 
quantities 2’ and 2” are negative when the 
suspension members and tires compress, since z 
is positive downwards. 

It will be useful in the following development 
of the theory to express motions of, or corre- 
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Fic. 2.’ Displacements at wheel “‘1.”’ 


sponding to, individual wheels in terms of net 
vertical translation (‘‘bounce’’), pitch, roll, and 
warping. Let us first consider motion of the 
ground, as defined by vertical motions of the 
wheel contact points. Referring to Fig. 3, let 
P,, P2, Ps, and P, be the contact points of the 
four wheels, which have risen by amounts 2), 22, 
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Z3, and z, above the plane of reference (undis- 
torted ground plane). Let NV, N’ be the midpoints 
of P:P, and P2P; and M, M’ be the midpoints 
of P,P: and P3P,, respectively. Roll is then 
defined as the slope of 14M’, pitch as the slope 
of NN’, and warp as half the relative inclination 
between P,P, and P2P3. We then have: 
Vertical translation of ground 


Zo = (21 +22+23+24)/4. (2) 
Roll o = (21+22—23—24)/2T, (3) 
where 7 = track. 
Pitch 09 = (2: —%2—23+24)/2W, (4) 
where W=wheelbase. 
Warp Wo= (21—22+23— 24) /2T, (S) 


¥» being 3 the warp angle. The roll, pitch, and 
warp angles are those of the distorted ground 
relative to the undistorted reference plane 
(angles being taken as equal to their tangents). 

Similarly, if zo’, Po’, 00’, wo’, are, respectively, 
the vertical translation, roll, pitch, and warp of 
the wheel center plane relative to the actual 
ground plane, these can be expressed in terms of 
tire deflections 2,’, 22’, 23’, zs’ by equations of the 
same form as (2), (3), (4), and (5). 

The net deformations of the suspension are 
similarly expressed in terms of 2,’’, 22’, 23’ 
these deformations being 29’’, o’’, 60’, Wo’ for 
vertical translations, roll, pitch, and 
respectively. 

The net motions of the suspended mass, which 
are denoted by 29’, Bo’, 00°", Wo” for vertical 
translation, roll, pitch, and warp, respectively, 


vr? vn? vrr 
9 9 


are similarly expressed in terms of 2;'’, 29’’, z 
z,"". The chassis is assumed to be perfectly 
rigid in torsion, and therefore Yo’ =0. 

The relation the motions of the 
ground, tires, suspension, and body have been 
expressed for each wheel in Eq. (1). Similar 
relationships apply to the net motions, and 
hence we have: 


o 
9 #4 +5 


warp, 


’ 


between 


Zo’ =Zo +21 +20”, 
,/"" = Pot+ho' +h”, 
6o”" = 69+ 0o'+ 89", 
Yo” a Voto +o” =0 


(6) 
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Fic. 3. Angular displacements. 


Loads in Suspension Members 


We have said that-in the most general type of 
suspension the load in the suspension member 
corresponding to any one wheel is affected by 
the displacements of every other wheel. Hence 
we can write for the load on wheel 1: 


Fy = — (A 421" + Byze!’ + C23" +Di24"’) (7) 


(loads being assumed proportional to deflections), 
where as previously stated 2;’’, 20’’, 23'’, 24’ 
measure the deflections of the four wheel centers 
relative to the sprung mass. 

The above equation can be written 


— Fy = F(2y" +29" +23" +2,'’) 
+ R(21" +22" —23" — 2,4" 
+P(21"’ — 20" —23" +24'") 
+H (2," —29''+23'’—2,4'"), (8) 
where 
A,=F+R+P+H, 
B,=F+R-—P-H, (9) 
C,= F—R-—P+H, 
D,=F-—R+P-H. 


Similar expressions, with the signs changed 
where relevant, apply to the other wheels. 
Thus for wheel 2, the coefficients in the expres- 
sion corresponding to (7) will be As, Bs, Cs, De, 
where 


A,=F+R-—P-H, ete. 


(Equation (8) and corresponding expressions are 
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always valid since the determinant 


a a 





‘4 4 —-1 1 
1 -1 -1 1 
‘4-1 1 —1 


is not zero.) 

Now the expressions in brackets in Eq.-: (8) 
are precisely the expressions given in the 
previous paragraph for the net deformation of 
the suspension in bounce, roll, pitch, and warp, 
respectively, and Eq. (8) can therefore be 
written: 


= Fy => 4F2,''+2RT," 
+2PW6)"+2HTYy"’. (10) 


From a consideration of this expression the 
coefficients F, R, P, H acquire a precise physical 
significance, since it is apparent that they 
represent the rigidity of the whole suspension in 
bounce, roll, pitch, and warp, respectively. 
These four rigidities are written £, \, u, v, where 


—t=—-—4F=dV/dzy", 

—\= —4RT?=dL/d4y", 11) 
—y=—4PW?=dM dy", 
—y=—4HT =2dN/dWo", 


where V, L, Mo, N are the net vertical load, 
rolling moment, pitching moment, and warping 
moment, respectively, acting on the whole 
suspension. 
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It may be explained for the sake of clarity the suspended mass (which is assumed torsionally 

that “warping” is a motion where the wheels, Tigid) except insofar as it sets up torsional 
: . rv . ” 

taken in rotation, go up, down, up, down. The Stesses. Nevertheless this Daten 2% has an 

; important effect on the motions of the wheels, 

warping of the ground has no direct effect on ' 

as will be shown later. 


Equations of Motion for Wheels 


The motion of the ground whether measured by 20, Wo, Po, 0, or by 21, Z2, 23, Z4 is a function of 
the shape of the ground surface and of the speed of translation u of the vehicle. The rear wheels 
are assumed to follow the same tracks as the front wheels, and hence: 


22(t)=2,(t+W/u), 
23(t)=2,(t+W/u), (12) 
or 

zo(t+W u) —29(t)= —-W 2(00(t) + 00(t+ W u)), 


D)(t+ W u) —Po(t) = —1 2(Wo(t) +Wo(t+ W. u)). ae) 


It follows that we can at will consider one or the other of the series of parameters 21, 22, 23, 24 
or 20, Wo, Bo, %, which may be considered as functions of time given in advance. The latter set will 
be chosen as it will lead to results which have a more direct physical interpretation. 

The positions of the axes of roll and pitch and the c.g. of the vehicle have important effects if the 
motion is curvilinear or if the vehicle is accelerated. In the case of rectilinear motion (which we 
are considering), the height of the c.g. above the roll and pitch axes gives rise to secondary effects 
which are: (a) an anti-stabilizing force; (b) a side force (due to roll) and forward force (due to pitch). 
The first effect is very small and may be neglected, while the second is neutralized by the adherence 
of the tires (for roll) and the freedom in forward motion (for pitch), and hence may be also neglected. 

The total motion of wheel 1 is given by: 


Zot+Zo +1/2(T(Po+0' +Yot+Yo') +W(80+60')), (14) 


and similarly (with signs changed where necessary) for other wheels. 
The suspension load on the same wheel is given by: 


1 1s/A m v 
F,= _ n't —( "4+ 0W"+—W0" ). (15) 
4 8\T W T 


Let K be the elastic rate of the tires. Hence the ground load on wheel 1 is 
Kz,’ = K(29' +4} T(®o0' +o’) + Wy}. (16) 
If m is the mass of the wheel, we can now write the equation of motion of wheel 1 as follows: 


ad 
pace tae +4 | T (Po +o’ +Wot+ Wo’) + W(0o+60’) } ] 
dt? 
1yé 1/X im y T W | : 
=- 2" +( y+" +40") -K(s'+ (20+ 00) +—60') . (17a) 
m\4 8\T W T 2 2 
Each wheel will have an equation of motion similar in form to (17a), but differing in the sign of 
the angular components. Let (17b), (17c), (17d), be the equations for wheels 2, 3, and 4. Then 
by adding all the four equations we shall eliminate all components but “‘bounce’’; by adding (17a), 
(17b), —(17c), —(17d) we shall eliminate all components but roll, and similarly, mutatis mutandis, 
for pitch and warp. We then obtain equations of motion in the various modes for the wheel system 
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as a whole. These equations will be: 


d? 1 
4— [Zo+20’ |=—[20’"’ —4Ka2y’ ], (18) 
dt? m 
d* TX 1 
4—[T(#)+ 4%’) J=— —o''—4KT®,’ |, (19) 
dt? mL 
d* 1fu 1 
4—[ W(0)+60’) ]=—| —60" —4K Woy |, (20) 
dt? mLW J 
q ifv 7 
4—[T(Wot+Wo') J=—| —Wo" —4K Tyo’ J. (21) 
dt? mLT 7 








Our final solution will require the partial elimination of the variables 2, zo’, ®o, Bo’, 0, 80’, Wo, Wo’. 
This can be done by making use of Eq. (6), which it will be recalled relates the above variables 
(which represent ground motion and tire compression) to variables representing motion of the 
suspension and the suspended mass. We then get: 


d? 1yé 
: (g9""’ — 29’) = — sto" —Kes""+K(eut20")| (22) 
dt? mL4 
d? if A 
cre (do'"’ — by’) =— - nite" Ko" + K bot"), (23) 
dt? mL4T? 
d? if pu 
—(00""’ — 0") =— cate" Kits" +K (00+60") |, (24) 
dt? mL 4W? 

dy." 17 » 

_-——_ =— cle" + Koo"). (25) 
dt? m4W? 





Equations of Motion for Complete System 


Let M be the suspended mass and let r?M and p*M be its moments of inertia in roll and pitch, 
respectively, 7 and p being the radii of gyration about the respective axes. The equations of motion 
of the suspended mass will then be: 





d?zo/"" 1 

—vanon th anced (26) 
dt? M 

d?@,/"" 1 

——=——_&,", (27) 
dt? rM 

d*6,/"" 1 

— a — $9”. (28) 

dt? p?>M 


We may substitute for d?z9’’’/d? in Eq. (22) from Eq. (26), giving an equation which we shall 
call (22a). Equation (22a) and Eq. (26) form the ‘‘system”’ S;: 





d?z,"’ 1 1 K K K 
+[e(—+— ) +—|w"—— ay" = —— (22a) 
| dt? M 4m m M m 
Si } 
d7zo'" ¢& 
| ——+— 29" =0, (26) 
| dt? M 
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This “‘system” completely defines the ‘‘bounce”’ (i.e., mean vertical translation) of the suspended 
mass and the suspension in terms of the mean vertical rise 2) of the ground surface. 

Similarly we can derive systems S2, S3, and S, which will define roll, pitch, and warp of the sus- 
pended mass and suspension (suspension only in the case of warp) in terms of the roll, pitch, and 
warp of the ground plane. Thus: 


d*h"" 1 1 K K K 
+{a( “) + ) “fh J»"- p,/" = _— Po. (23a) 
| dt? 4mT? #r°M m m m 
So 
| d*@,/"" r 
+ ,'’=0, (27) 
dls r°M 
d*6,"’ 1 1 K K K 
+| + )+ Jor 0," = = Ao, (24a) 
| d# 4mW? p?M M m m 
3 Prva . 
| 7s 6,’’=0, (28) 
dt? p> M 
dp," v K K 
S4| ; +| — y= ——Ho. (25) 
| dt? 4mW2 m m 


S, is merely a repetition of Eq. (25). With the exception of this last system, the above are all 
of the same form and may be completely solved if the independent variable (which represents 
“‘motion”’ of the ground under the vehicle) is known as a function of time. 


STUDY OF MOTION 
General Remarks 


It will be shown below that the system which constitutes a four-wheel vehicle has seven natural 
modes of oscillation, each with its own frequency; two modes each correspond to bounce, roll, and 
pitch, and one to warp. The number of different ‘modes can be deduced a priori from the fact that 
the system has seven degrees of freedom, viz., the four compressions of the tires and the bounce, 
roll, and pitch of the body. 

To analyze the motion we must solve the equation systems S,; to S;. The general solution is the 
sum of the general solution with the right-hand sides zero and a particular solution with particular 
values of the right-hand sides. The general solution will give us the natural frequencies of the system, 
and the complete solution will give the total motion when passing over a bump of known shape. 
A representative shape of ground bump will be selected for this purpose, a procedure which is justified 
by the fact that the actual shape chosen has no influence on the general conclusions to be drawn 
from the results of the analysis. 

We shall begin by studying the system Sj, i.e., the motion in the ‘“‘bounce’’ mode. 


General Solution of Equations of Motion 
The general solution of system S, (Eq. (24), and Eq. (22a) with R.H.S.=zero) will be of the type 
p2o"’ +qz0'"’ =A sin (wt+e), 


and p and q are coefficients to be determined. Hence 


d?zo"’ d*Z 


p- —-+-¢——_ = —w?A (p29 +920’). 
dt* dt? 
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Substituting this in system S, to eliminate d*z9’’/dé and d?z9’"’/dé’, 


1 1 K K E 
pl (— +) +— od" — pa” 92!" —(pe" +40") =0. 
M 4m m m M 


From which we may equate coefficients of 29’ and 29’”’ to zero, getting: 


1 1 K E 
| :(—+—) +—|+¢——w'p =0, ee 
M 4m m M 
K 
—p——w'g=0. (27) 
m 
Hence we get 
1 1 K g 
M 4m m M 
=0, (27a) 
—_ — —w? 
m | 


which gives us two roots for w. Let these roots be ai, 6:. From Eq. (27) we have 


p w*m 


q K 


The two values for the general solution of system S, are, therefore 


m 

Zo’ —ay— 2p = 1, cos (a,;t+a,), (28) 
K 
m 

alll dina, cos (8;t+)). (28a) 


Particular Solution and Effect of Ground Bump 


Let us consider the effect of one wheel passing over a bump of maximum height A,,, this maximum 
being reached at a distance x from the start of the bump, the shape being represented by an expo- 
nential law, viz.: 


K 
—Z9 = Ti? exp —rl, (29) 
m 
where 
1 u 
r=-- (u=forward speed) 
ex 
and 
m y ™ 
4T—=u—, 
K x 


zy) being the mean rise of the ground, i.e., one-quarter of the height of the bump. 
If we substitute this in the right-hand side of Eq. (22a), we shall be able to find the particular 
solution of the system S;. It can be shown that the right-hand side of this solution will be the same 
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as for the general solution, and that the answer will be: 


a\"*m r 
Zo’ — Zo = (t+7,’) exp —ft, (30) 
K a;-+r* 
B\*m r 
Z,/"’ — Zy'= (t+7,") exp —rt, (30a) 
K Bi°+r* 
where 
2r s 2r 
wT, = ’ ) = —, 
ai?+r* Bi+r 


Motion of System 


The complete solutions of the system S; are the sums of (28) and (30), and (28a) and (30a). 
Separating the variables 29’ and 29’, we get the motions of the suspended mass and the suspension 


members relative to the suspended mass. The motion of the suspended mass is given by: 





1 
zy’ = | ie cos (a;f+a;) —Jip1" cos (B\t+;)} 
pi’ — pi’ , ” 
| pi P ’ pi P ” \] 
+exp—rt . (t+ 7, )-— (t+7, ) (31) 
las?+r? Bi2+r? 
and that of the wheels in relation to the mass by: 
1 r r \] 
Zo = as | 1: c0s (8,t+6,) —I; cos (a;t+a,)} +exp—rt: —(t+7,"")— - -(t-+71') (32) 
Pi —Pi lg.2+r? ay*+r* 
where 
m m 
pi” 


=a—, pi =p. 
K K 


The constants J;, Ji, a1, 6;, can be found from a knowledge of the initial conditions. If the system 
is assumed to be at rest at the instant of hitting the bump, then: 
a;"?—r* r 


tan a,;= —, [,=—————__,,_ et. (33) 
2air a;(ay*+r?) 


Motion in Roll, Pitch, and Warp 


The motion in roll is determined in exactly similar fashion as that for bounce, by solving the 
equations of system S:2, which is of the same form as system 5;. The oscillations in roll have angular 
velocities a2 and 82 (angular velocity = 27 X frequency) which are the roots of the equation 


» K 2» ~~ X 


+—t+ —w* — 
4m7T?> m PrM r°>M 
== (). 
K 
_ —w?" 
m 


The same ground bump which we have considered above (if occurring under wheel 1) is equivalent 
to a “rolling” of the ground plane given by 
K 2 


$)=—Tt exp—rt. (34) 
m 7 
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The complete solution of system S2 will then give the rolling motions of the suspended mass 
(#,'"") and of the wheels in relation to the mass (# 9’) in a form identical with that of Eqs. (31) 
and (32). The suffix 2 will be used for the constants involved in roll. 

The motion for pitch is obtained in similar fashion by solving equation system S3, the natural 
angular velocities in this case being a3 and 83. 

The solution of Eq. (25) or S, will give the wheel motion in warp (the suspended mass cannot 
warp), the general solution being: 


Wo’ = —I, cos (agt+as), (35) 


y Ky} 
a= (— +) ° 
4mT? m 


Net Motion 


where 





We have now analyzed the seven distinct vibrations which are excited in the wheels and suspension 
by the passage of one wheel over a given obstacle. Except for the seventh (warp), these vibrations 
are all transmitted in various degrees to the suspended mass. The vibrations contain a sinusoidal 
component and a non-periodic component, the latter representing the shape of the obstacle, which 
is reproduced in the motions of the suspension and the sprung mass with some reduction of amplitude 
and change in phase. 

These motions may be grouped as follows, for the suspended mass: the three principal vibrations 
correspond to angular velocities 81, 82, 83, and to bounce, roll, and pitch, respectively. The three 
vibrations corresponding to angular velocities a1, a2, a3; (and also to bounce, roll, and pitch), may 
be regarded as terms due to coupling of the suspended mass with the wheels. The only effect of 
the warping vibration (a4) on the suspended mass is to induce torsional stresses in the chassis. 

Similarly for the wheels the terms in a1, a2, a3, as, may be regarded as the principal vibrations 
and the terms in 8), 82, 83 as coupling effects with the suspended mass. 

The net motion of wheel 1 (relative to the body) is given by: 








T W T 
3)’ = 29 +—P 9" +— 60" +—Wo"”, (36) 
: 2 2 2 
i.e., 
1 T 
2,;'= Re >t Jicos (B:t+6,) —I, cos (a;t+a,)}+— enters | J2cos (B2t+b2) — I2cos (aet+aze)} 
Pi —Pi 2(po’’ — po’) 
W T 
+————| J; cos (83t-+63) —I3 cos (a3t+a3) } ——I,4 cos (agt+a,4) 
2(ps'’ — ps’) 2 
rexp—ri | 1 1 1 
sai “| a OS all ” P a tea alee a rene 26 tin ae 
2 (Bi°+r*)(pi —pi) (B2°+9r*)(p2'—p2) (Bs’-+1r*) (ps — ps) 
1 1 1 1 








(axt+r2)(pi""—pr') (aat+r*)(pe"— ps’) (as*+r°)(ps" pa’) (aat-+r2) 


,? ”? ” I. 
T1 T2 T3 T1 








+ -— meee = aaa meaner 
(Bi?+r?)(pi’— pr’) (B2®+1r?)(p2"”— po’) (Bs?-+1?)(ps’"’— ps’) (ar? +r?) (pi — py’) 





U LU / 





T2 T3 T4 
(as?-+r*)(p2""— pa!) (as*+r*)(ps"— ps!) aa?+r2d 





(37) 
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It should be noted that the elimination of the parameter », the rigidity in warp (i.e., making this 
equal zero) will not modify the form of the preceding results; its only effect would be to eliminate 
torsional stresses on the chassis and to reduce ay. 

Integration of the equations of motion can also be performed for variations of ground shape of 
sinusoidal form, or form which may be expressed in terms of sinusoidal functions. 


DISCUSSION OF RESULTS 
Orthodox Suspensions 


In the case of orthodox suspensions the relative 
rigidities of the suspension are interrelated by 
laws which are affected only by the proportions 
of the chassis (to a small extent). This means 
that there is a relation between the four coeffi- 
cients F, R, P, H or the four quantities A,, Bi, 
C,, D,; and their homologues. 

For example in the case of an orthodox 
suspension with axles there is no interaction 
between wheels on the same side of the vehicle, 
and therefore: 


B,=0, A.=0, A;=0, B,=0, 


D;=0, C,=0, 


which gives F+ R—P—H=0, and F—R—P+H 
=0, iec., R=H and F=P, whence \=v and 
&é=p/W?. 

We shall show that for suspensions of this 
type the only coupling between wheels which 
are not on the same side is through the suspended 
mass. Considering the bounce mode, the angular 
velocities a; and 6; in this mode are given by 
Eq. (27a). If now we make M (the suspended 
mass) infinite in this equation, the root a 
becomes a,’ = (£/4m+K/m)'. Similarly for pitch, 
the root a3 becomes a3’ = (u/mW?+K/m)}. Hence 
a, and as are equal, but the signs of the vibra- 
tions corresponding 


C,=0, D,=0, 


to a; will be opposite for 
wheels on opposite axes. Thus if one of the 
front wheels passes over a bump the bounce 
and pitch oscillations in the rear wheels will be 
of the same magnitude but opposite sign (i.e., 
out of phase by 180 degrees), and the net 
oscillation in the rear wheels due to pitch and 
bounce will be zero. Hence the only coupling 
must be through the suspended mass, since, 
when this mass is made infinite, the coupling 
disappears. 

In the case of a symmetric car with four 
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(Note: Coefficients J; and J; have one dimension in length. Coefficients J2, Jz, etc., are non-dimen- 





independent wheels, the only coupling between 
any wheels is through the intermediary of the 
suspended mass. In this case F=P=H=R. 


Causes of Discomfort in Orthodox Suspensions 


Due to coupling effects the seven frequencies 
of vibration are always distinct. However, the 
angular velocities a1, a, a3, as, are in fact very 
close in value, and this is the reason for the 
‘‘beat’”” phenomena which are the main cause of 
discomfort in cars with orthodox suspensions if 
there is inadequate damping. 

“Beats” occur if there are two superimposed 
oscillations, the frequencies of which are fairly 
close. We shall now proceed to study this in 
more detail, and to determine the degree of 
closeness between some of the frequencies and 
the factors which influence it, for the case of an 
orthodox independent wheel suspension. 

In this case we have F=R, and therefore 
A= T°. 

The factors 1/M+1/4m and 1/Mr?+1/4mT? 
occur in the expressions for bounce and roll 
frequencies. We may write the latter expression 


as 
1 /T? 1 1 
( 2 * , ) , 
T?\r? M 4m 


The radius of gyration in roll 7 being normally 
less than the track 7, we may write 


T? 1 1 1 1 
+= +, 
r> M 4m M 4m 


where ‘“‘x’’ is greater than unity. 
Now the vibration frequencies in bounce and 


roll are given by: 
1 1 K g 
(—+—)+—-e, _ 
M 4m m m 


m 






— 
’ Ww 
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and 


|X 1 1 K A | 
u(— +) +- —«*, —| 
T? \M 4m m r? M | 
| 


| K 
| ences o--— 
| m | 


the last of which (roll) may be written: (since 
\/T?= 8) 


| 1 1 K u’t 
eu( + -) += —w, - 
| M 4m m M 


'=0, (38) 
K 
_--—-, —w" 
m 
where 
T? 
u’=—, (u’>1), 
r- 
i.e., 
1 1 K EK 
ote tu —-+—)-+--|+u'—=0, (39) 
M 4m m Mm 


This equation will give the roll frequencies by 
inserting the proper values for u and uw’, and 
also the bounce frequencies by making u and w’ 
equal to unity. 

Increasing the rigidity in roll independently 
of that in bounce is equivalent to increasing u 
and u’ in Eq. (39), i.e., to increasing both the 
sum and product of the frequencies in roll as 
compared with bounce. The effect therefore 
would be to. “decouple” the bounce and roll 
oscillations and to reduce beats. 


In an axle type of suspension \<£/T and 
hence u and wu’ tend to approach unity. Hence 
the frequencies for roll and bounce are very 
close together. This is one of the main reasons 
for the advantage of independent suspensions 
over those of the axle type. 


“Interaction” Suspensions 


Our study has covered the most general type 
of suspension of which: the axle type and the 
independent wheel type are but particular 
instances. It is apparent from what has been 
said above that the independent wheel sus- 
pension has no fundamental advantages from a 
general standpoint, although it may be preferable 
to the axle type. 

The fundamental solution is a suspension in 
which the four quantities F, R, P, H, or in other 
words the rigidities in bounce, roll, pitch, and 
warp, can be selected entirely independently of 
each other, thus allowing the choice of the best 
compromise between stability and comfort (a 
step in that direction is the introduction of roll 
stabilizers). In particular, it is possible to 
eliminate rigidity in warping and thus suppress 
torsion loads on the chassis. 

Such suspensions have actually been devised 
and tested by the author. They have proved to 
be very much superior to any other type, both 
from the point of view of comfort and road 
holding; in fact, it has been possible in some 
cases to dispense with shock absorbers without 
detrimental effect on comfort. 
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Calculation of Axially Symmetric Fields* 





SIDNEY BERTRAM 
The Ohio State University, Columbus, Ohio 


(Received January 9, 1942) 


A practical method is developed for calculating the 
distribution of scalar fields in axially symmetric systems. 
To obtain a complete solution a knowledge of the potential 
distribution along a cylinder of constant radius is required; 
in many cases the distribution along the cylinder can be 
estimated with sufficient accuracy by an inspection of the 
electrode configuration. The solution is first obtained in 
the form of an integral. This is evaluated for the class of 
boundary potential functions of the form 
V.(1, 2) = s*, s>0 

0,2<0. 
In this expression, V,(1, 2) is the potential along a cylinder 
of unit radius, z is measured along the axis of the cylinder, 


HE increasing number of devices in which 
scalar fields with symmetry about a central 
axis occurs makes it desirable to have a simple 
method of calculating the field in such cases. 
In particular, such devices as kinescopes, image 
converters, etc., require a knowledge of the 
potential fields before a calculation of the elec- 
tron trajectories can be made. In the past 
electrolytic plotting tanks have been used to 
obtain equipotential maps for particular elec- 
trode systems and a graphical method is then 
used to obtain the electron trajectories. The 
accuracy of such calculations is limited by the 
refinement of the apparatus used to obtain the 
potential plot. This paper gives a_ practical 
method for calculating the field in such cases. 
Scalar potential fields must satisfy Laplace's 
equation; for axially symmetric fields it is con- 
venient to consider this equation expressed in 
cylindrical coordinates: 


aV 10aV aV 
v2V=—+ +—=0, (1) 


Or? ror dz" 


where V= V(r, 2) is the potential at a radius r 
from the axis of symmetry and at a distance z 
from a reference plane perpendicular to the axis. 
Several general solutions of this equation have 


* Taken from a thesis of the same title presented in 
partial fulfillment of the requirements for the degree of 
Master of Science in Electrical Engineering (1941) at 
The Ohio State University. 
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and m is a positive integer (zero included). Numerical 
values are given in tabular form for the cases n=0, 1, 2 
and 3. It is then shown how the solution for the above 
boundary potential functions can be used to obtain the 
solution for an arbitrary distribution of potential along a 
cylinder of constant radius. In particular, the potential 
distribution of a simple electron lens consisting of two 
coaxial cylindrical electrodes of equal diameters is expressed 
in terms of the tabulated functions. The calculated equi- 
potential lines for one such case are compared graphically 
with a plot made in an electrolytic plotting tank. The 
comparison between the calculated and experimental 
results is quite favorable. 


been given;! the particular solution used here is 
the Fourier integral 
1 a 
V(r, s) = -f a(k)Io(krje**dk, (2) 
2n J_, 
where Jo(kr)=Jo(ikr) is the Bessel function of 
order zero with imaginary argument and where 
a(k) is a function of the arbitrary parameter k; 
a(k) must be chosen so that V(r, 2) satisfies the 
boundary conditions. 
If the potential at some radius R is known asa 
function of z then a(k) can be found by means of 
the Fourier transform 


a(k)Io(RR) -{ V(R, z)e~***dz (3) 
or, if for convenience the radius at which the 
potential is known is taken as the unit of linear 
measure, then R=1 and (3) becomes 


x 


a(k)Tu(k)= f V(1, s)e~ “dz. (4) 


—@® 


From (4) a(k) can be found explicitly in terms 
of the boundary potential V(1, z). Substitution 
of a(k) thus found in (2) yields 


; 1 as 24 Io(kr) 
Vir, 2) = -f f Vii, ¢)——e*“-Ddkd¢, (5) 
2 —o Yo To(k) 


' See, for example: L. M. Myers, Electron Optics, Theo- 
retical and Practical (D. Van Nostrand Company, New 
York, New York 1939), pp. 87-91. 
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where the variable of integration in (4) is re- 
placed by ¢ to avoid confusion. Since the order 
of the integration in (5) is immaterial,? the po- 
tential can be expressed 


® 


V(r, =f V1, U(r, de, (6) 
where ss 
1 ” Io(kr) 
U(r, §) = -f —_—eittdk (7) 
and = 


Since (7) is independent of the boundary 
potential V(1,z) it can be evaluated and the 
result used in (6) to obtain the potential distribu- 
tion for the general case. To evaluate (7) the 
fraction Jo(kr)/Io(k) is first expanded into a 
series of partial fractions,‘ thus 


Io(kr) To(tyr) So(ur) 
—= 5 ——_____ =), (8) 
Io(k) = * (R—tp)Io'(ip)  * iTs(u)(k —in) 





where the summation is taken over all roots yp of 
the equation Jo(u) =0. Then 


1 *  Jo(prje* 
Ue=—E J 5O dk. (9) 
2ni * Js(u)(k—in) 





Evaluating (9) by a contour integration there 
results: 

* The order of the integration in (5) is immaterial if the 
integral (4) converges. This is not true for boundary po- 
tentials V(1, 2) as defined in (14). The condition is satisfied, 
however, for functions of the form 


(0, 2s>A 
V,*(1,2)= 42", O<2<A 
(0, 2<0 


This would require that the upper limit of integration in 
(15) be replaced by A. Since (15) converges uniformly as 
A-—»« the results here obtained by changing the order of 
integration in (5) are valid. 

’ For the calculation of paraxial electron trajectories the 
axial potential distribution is sufficient. A good approxima- 
tion for this purpose is 


U(0, ¢) = 5 sech? wé, 
where w= 1.32 giving 


V(0, z) ~SS_ V(1, £) sech? wide. 


For a discussion of this equation see S. Bertram, ‘‘Deter- 
mination of the axial potential distribution in axially 
symmetric fields,’’ Proc. I.R.E. 28, 418-420 (1940). 

*B.O. Peirce, A Short Table of Integrals (Ginn and Com- 
pany, 1929), pp. 13-14. 
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A. for positive pu 


1 © pikt et £>0 
— f inal -dk= | , (10) 
0 


2ni J _. R—ip 


B. for negative p 


1 - @ 0, &>0 
—{ t= (11) 
2nid_. k—ip e® E<0. 
Since there is a negative root numerically equal 


to every positive root of Jo(x)=0 and vice 
versa (9) can finally be written 


J o(ur) 
U(r, t) = etl, (12) 
+h J(u) 


In terms of the series (12) Eq. (6) becomes 


J o(ur) 
Ve,9<E — Fay f Vee wena. a9) 
+u Ji Le D 


It is now expedient to apply (13) to the simple 
class of boundary functions of the form 


0, 2<0 


V,(1, 2) = (14) 





2", 2>0, 


where n is a positive integer or zero. This yields 
Vl. )= Ef penenas (15) 


For z>0 this becomes 


Jo(ur) 
V(r, 2) => - Af gre #e-Ddg 


z>0 +u J(u) 


+f renee 


Jo(ur) | , 
=> ‘ui Hf err(- —----—- 


+u J;(u) Mm pn 


ce” 2 z 
+n(n—1)——+-:: 
p? 


—¢* nor” 1 
er] e-ot(— ee eens 
rm a 


tient ae . >) | (16) 
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TABLE I. 2 o(ur) gn-2 
V,(r, 2) =2"+ >> ————j n(n—1)—— 
z>0 <0 2>0 +u Ji(u) pn 
Volr, 2) 1— > Jo(ur), Bz >> Jour) ous ili 
ru mS i() +m BS (us) +n(n—1)(n—2)(n—3)—— 
A 5 
, as y Jo(ur) Me y Jo(ur) " KM 
Vilr, 2) Se Adal Say 
- . . lp—uz 
2Jo(ur) 2Jo(ur) nN € o 
V2(r, 2) +> — [2-e#] 3" eee}, (17) 
tu mS (pe) iu mJ 1 (p) 2(—p)"*! 
Valr, 2) 33 > sam ke 4° | y 6Jolur) ous I Or 2 <0 r 
tu me? JS () 7 ru pid (pm) r J o(ur) 
f V(r, 2) => of cree 
2<0 +m Ji(p) . 
. , ‘ \a lpuz 
One obtains from Vo(r, z) at z=0 the relation af oo satel (18) 
+ ys ( ) n+l 
2Jo(ur) — 
te pJi(ps) Table I lists the first four of these functions (nm =0, 
1, 2, 3). The series are all convergent; the more 
using this (16) reduces to important values are tabulated in Tables II-V. 
Tapee I. > 28) ¢-us 
tu pd (mu) 
Nr 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
0.00 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 
0.05 0.4669 0.4666 0.4656 0.4639 0.4610 0.4566 0.4495 0.4372 0.4126 0.3441 
0.10 0.4341 0.4335 0.4316 0.4281 0.4226 0.4140 0.4005 0.3778 0.3347 0.2360 
0.15 0.4019 0.4010 0.3983 0.3933 0.3854 0.3732 0.3545 0.3240 0.2706 0.1697 
0.20 0.3706 0.3695 0.3660 0.3596 0.3497 0.3348 0.3121 0.2769 0.2201 0.1279 
0.25 0.3404 0.3391 0.3349 0.3276 0.3161 0.2991 0.2740 0.2367 0.1808 0.1000 
0.30 0.3116 0.3101 0.3054 0.2973 0.2847 0.2664 0.2401 0.2027 0.1501 0.0803 
0.35 0.2841 0.2825 0.2776 0.2689 0.2557 0.2367 0.2103 0.1740 0.1259 0.0658 
0.40 0.2583 0.2567 0.2515 0.2425 0.2290 0.2100 0.1842 0.1499 0.1065 0.0548 
0.45 0.2342 0.2325 0.2273 0.2182 0.2047 0.1861 0.1614 0.1296 0.0907 0.0461 
0.50 0.2118 0.2101 0.2049 0.1959 0.1827 0.1648 0.1416 0.1124 0.0778 0.0392 
0.55 0.1910 0.1894 0.1843 0.1755 0.1629 0.1459 0.1243 0.0978 0.0671 0.0336 
0.60 0.1719 0.1703 0.1654 0.1571 0.1451 0.1292 0.1092 0.0853 0.0581 0.0289 
0.65 0.1545 0.1529 0.1482 0.1404 0.1291 0.1143 0.0961 0.0746 0.0504 0.0250 
0.70 0.1385 0.1371 0.1327 0.1253 0.1148 0.1012 0.0846 0.0653 0.0440 0.0217 
0.75 0.1241 0.1227 0.1186 0.1117 0.1020 0.0896 0.0745 0.0573 0.0384 0.0189 
0.80 0.1110 0.1097 0.1059 0.0995 0.0906 0.0793 0.0657 0.0503 0.0337 0.0165 
0.85 0.0991 0.0979 0.0944 0.0886 0.0805 0.0702 0.0580 0.0443 0.0295 0.0145 
0.90 0.0885 0.0874 0.0841 0.0788 0.0714 0.0622 0.0512 0.0390 0.0260 0.0127 
0.95 0.0789 0.0779 0.0749 0.0701 0.0634 0.0551 0.0453 0.0344 0.0228 0.0112 
1.00 0.0703 0.0694 0.0667 0.0623 0.0563 0.0488 0.0400 0.0303 0.0201 0.0098 
1.05 0.0626 0.0618 0.0593 0.0554 0.0499 0.0432 0.0354 0.0268 0.0177 0.0087 
1.10 0.0557 0.0550 0.0527 0.0492 0.0443 0.0383 0.0313 0.0237 0.0157 0.0076 
1.15 0.0495 0.0489 0.0469 0.0437 0.0393 0.0339 0.0277 0.0209 0.0138 0.0067 
1.20 0.0440 0.0434 0.0417 0.0388 0.0349 0.0301 0.0245 0.0185 0.0122 0.0059 
1.25 0.0391 0.0386 0.0370 0.0344 0.0309 0.0266 0.0217 0.0164 0.0108 0.0053 
1.30 0.0347 0.0343 0.0329 0.0306 0.0274 0.0236 0.0192 0.0145 0.0096 0.0047 
1.35 0.0309 0.0304 0.0292 0.0271 0.0243 0.0209 0.0170 0.0128 0.0084 0.0041 
1.40 0.0274 0.0270 0.0259 0.0241 0.0216 0.0186 0.0151 0.0113 0.0075 0.0036 
1.45 0.0243 0.0240 0.0230 0.0213 0.0191 0.0164 0.0134 0.0101 0.0066 0.0032 
1.50 0.0216 0.0213 0.0204 0.0189 0.0170 0.0146 0.0119 0.0089 0.0059 0.0029 
1.55 0.0192 0.0189 0.0181 0.0168 0.0150 0.0129 0.0105 0.0079 0.0052 0.0025 
1.60 0.0170 0.0168 0.0160 0.0149 0.0133 0.0115 0.0093 0.0070 0.0046 0.0022 
1.65 0.0151 0.0149 0.0142 0.0132 0.0118 0.0101 0.0082 0.0062 0.0041 0.0020 
1.70 0.0134 0.0132 0.0126 0.0117 0.0105 0.0090 0.0073 0.0055 0.0036 0.0017 
1.75 0.0119 0.0117 0.0112 0.0104 0.0093 0.0080 0.0065 0.0049 0.0032 0.0015 
>1.75 0.800a 0.788a 0.754a 0.700a 0.625a 0.537a 0.434a 0.3260 . 0.214a 0.104 
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+u wd (mu) 
2\r 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
0.00 0.2674 0.2658 0.2607 0.2522 0.2398 0.2232 0.2016 0.1739 0.1377 0.0883 
0.05 0.2432 0.2416 0.2366 0.2281 0.2158 0.1993 0.1779 0.1505 0.1150 0.0674 
0.10 0.2207 0.2191 0.2141 0.2058 0.1937 0.1775 0.1567 0.1301 0.0964 0.0531 
0.15 0.1998 0.1982 0.1934 0.1852 0.1735 0.1579 0.1378 0.1126 0.0813 0.0431 
0.20 0.1805 0.1790 0.1743 0.1664 0.1551 0.1402 0.1212 0.0976 0.0691 0.0358 
0.25 0.1627 0.1613 0.1568 0.1492 0.1385 0.1243 0.1065 0.0848 0.0591 0.0301 
0.30 0.1465 0.1450 0.1408 0.1336 0.1235 0.1102 0.0937 0.0738 0.0509 0.0256 
0.35 0.1316 0.1302 0.1262 0.1195 0.1100 0.0976 0.0824 0.0644 0.0440 0.0220 
0.40 0.1180 0.1168 0.1130 0.1067 0.0979 0.0865 0.0726 0.0564 0.0382 0.0190 
0.45 0.1057 0.1045 0.1011 0.0952 0.0870 0.0766 0.0640 0.0494 0.0333 0.0165 
0.50 0.0946 0.0935 0.0903 0.0849 0.0774 0.0678 0.0564 0.0433 0.0291 0.0143 
0.55 0.0845 0.0835 0.0805 0.0756 0.0687 0.0601 0.0498 0.0381 0.0255 0.0125 
0.60 0.0754 0.0745 0.0718 0.0673 0.0611 0.0532 0.0439 0.0335 0.0224 0.0110 
0.65 0.0673 0.0664 0.0640 0.0598 0.0542 0.0471 0.0388 0.0295 0.0196 0.0096 
0.70 0.0600 0.0592 0.0569 0.0532 0.0481 0.0417 0.0343 0.0261 0.0173 0.0085 
0.75 0.0534 0.0527 0.0507 0.0473 0.0427 0.0370 0.0303 0.0230 0.0152 0.0075 
0.80 0.0475 0.0469 0.0451 0.0420 0.0379 0.0328 0.0268 0.0203 0.0134 0.0066 
0.85 0.0423 0.0417 0.0401 0.0373 0.0336 0.0290 0.0237 0.0179 0.0119 0.0058 
0.90 0.0376 0.0371 0.0356 0.0332 0.0298 0.0257 0.0210 0.0159 0.0105 0.0051 
0.95 0.0334 0.0330 0.0316 0.0294 0.0265 0.0228 0.0186 0.0140 0.0093 0.0045 
1.00 0.0297 0.0293 0.0281 0.0261 0.0235 0.0202 0.0165 0.0124 0.0082 0.0040 
1.05 0.0264 0.0260 0.0249 0.0232 0.0208 0.0179 0.0146 0.0110 0.0072 0.0035 
1.10 0.0234 0.0231 0.0221 0.0206 0.0185 0.0159 0.0129 0.0097 0.0064 0.0031 
1.15 0.0208 0.0205 0.0196 0.0183 0.0164 0.0141 0.0115 0.0086 0.0057 0.0028 
1.20 0.0185 0.0182 0.0174 0.0162 0.0145 0.0125 0.0101 0.0076 0.0050 0.0024 
1.25 0.0164 0.0162 0.0155 0.0144 0.0129 0.0111 0.0090 0.0068 0.0044 0.0022 
1.30 0.0145 0.0143 0.0137 0.0127 0.0114 0.0098 0.0080 0.0060 0.0039 0.0019 
1.35 0.0129 0.0127 0.0122 0.0113 0.0101 0.0087 0.0071 0.0053 0.0035 0.0017 
1.40 0.0114 0.0113 0.0108 0.0100 0.0090 0.0077 0.0063 0.0047 0.0031 0.0015 
1.45 0.0102 0.0100 0.0096 0.0089 0.0080 0.0068 0.0055 0.0042 0.0027 0.0013 
1.50 0.0090 0.0089 0.0085 0.0079 0.0071 0.0060 0.0049 0.0037 0.0024 0.0012 
>1.50 0:332a 0.327a 0.313 0.291a 0.260a 0.223a 0.180a 0.135a@ 0.089a 0.043a 
a= e242 
om >) J o(ur) ie 
TABLE IV. Saw e 
rw 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
0.00 0.1250 0.1237 0.1200 0.1137 0.1050 0.0937 0.0800 0.0637 0.0450 0.0237 
0.05 0.1122 0.1111 0.1076 0.1019 0.0936 0.0832 0.0705 0.0556 0.0387 0.0199 
0.10 0.1006 0.0996 0.0963 0.0910 0.0834 0.0738 0.0622 0.0487 0.0334 0.0169 
0.15 0.0901 0.0891 0.0861 0.0812 0.0742 0.0654 0.0548 0.0426 0.0290 0.0145 
0.20 0.0806 0.0797 0.0769 0.0724 0.0660 0.0580 0.0484 0.0373 0.0252 0.0125 
0.25 0.0721 0.0712 0.0687 0.0645 0.0587 0.0514 0.0427 0.0328 0.0221 0.0109 
0.30 0.0644 0.0636 0.0613 0.0574 0.0521 0.0455 0.0377 0.0288 0.0193 0.0095 
0.35 0.0574 0.0567 0.0546 0.0511 0.0463 0.0403 0.0333 0.0254 0.0169 0.0083 
0.40 0.0512 0.0505 0.0486 0.0454 0.0411 0.0357 0.0294 0.0224 0.0149 0.0073 
0.45 0.0456 0.0450 0.0433 0.0404 0.0365 0.0317 0.0260 0.0197 0.0131 0.0064 
0.50 0.0406 0.0400 0.0385 0.0359 0.0324 0.0280 0.0230 0.0174 0.0115 0.0056 
0.55 0.0361 0.0356 0.0342 0.0319 0.0287 0.0248 0.0203 0.0154 0.0102 0.0050 
0.60 0.0321 0.0317 0.0304 0.0283 0.0255 0.0220 0.0180 0.0136 0.0090 0.0044 
0.65 0.0286 0.0282 0.0270 0.0252 0.0226 0.0195 0.0159 0.0120 0.0079 0.0039 
0.70 0.0254 0.0250 0.0240 0.0223 0.0201 0.0173 0.0141 0.0106 0.0070 0.0034 
0.75 0.0225 0.0222 0.0213 0.0198 0.0178 0.0153 0.0125 0.0094 0.0062 0.0030 
0.80 0.0200 0.0197 0.0189 0.0176 0.0158 0.0136 0.0111 0.0083 0.0055 0.0027 
0.85 0.0178 0.0175 0.0168 0.0156 0.0140 0.0120 0.0098 0.0074 0.0049 0.0024 
0.90 0.0158 0.0156 0.0149 0.0139 0.0124 0.0107 0.0087 0.0065 0.0043 0.0021 
0.95 0.0140 0.0138 0.0132 0.0123 0.0110 0.0095 0.0077 0.0058 0.0038 0.0019 
1.00 0.0124 0.0123 0.0117 0.0109 0.0098 0.0084 0.0068 0.0051 0.0034 0.0016 
>1.00 0.1382 0.136a 0.130a 0.121a 0.108a 0.093a 0.075a 0.056a 0.037a 0.018a 
a= e242 
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Whites te Oe ma 
raB_e V. & Sin) 
Nr 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
0.00 0.0550 0.0543 0.0523 0.0490 0.0445 0.0389 0.0323 0.0248 0.0168 0.0083 
0.05 0.0490 0.0484 0.0466 0.0436 0.0395 0.0345 0.0285 0.0218 0.0146 0.0072 
0.10 0.0437 0.0431 0.0415 0.0388 0.0351 0.0305 0.0252 0.0192 0.0128 0.0063 
0.15 0.0389 0.0384 0.0369 0.0345 0.0312 0.0271 0.0223 0.0169 0.0113 0.0055 
0.20 0.0347 0.0342 0.0329 0.0307 0.0277 0.0240 0.0197 0.0149 0.0099 0.0049 
0.25 0.0308 0.0304 0.0292 0.0272 0.0246 0.0212 0.0174 0.0132 0.0087 0.0043 
0.30 0.0274 0.0271 0.0260 0.0242 0.0218 0.0188 0.0154 0.0117 0.0077 0.0038 
0.35 0.0244 0.0241 0.0231 0.0215 0.0193 0.0167 0.0136 0.0103 0.0068 0.0033 
0.40 0.0217 0.0214 0.0205 0.0191 0.0172 0.0140 0.0121 0.0091 0.0060 0.0029 
0.45 0.0193 0.0190 0.0182 0.0169 0.0152 0.0131 0.0107 0.0081 0.0053 0.0026 
0.50 0.0171 0.0169 0.0162 0.0150 0.0135 0.0116 0.0095 0.0071 0.0047 0.0023 
0.55 0.0152 0.0150 0.0144 0.0133 0.0120 0.0103 0.0084 0.0063 0.0042 0.0020 
0.60 0.0135 0.0133 0.0127 0.0118 0.0106 0.0091 0.0074 0.0056 0.0037 0.0018 
0.65 0.0120 0.0118 0.0113 0.0105 0.0094 0.0081 0.0066 0.0049 0.0033 0.0016 
0.70 0.0106 0.0105 0.0100 0.0093 0.0083 0.0072 0.0058 0.0044 0.0029 0.0014 
0.75 0.0094 0.0093 0.0089 0.0083 0.0074 0.0064 0.0052 0.0039 0.0026 0.0012 
0.80 0.0084 0.0082 0.0079 0.0073 0.0066 0.0056 0.0046 0.0034 0.0023 0.0011 
0.85 0.0074 0.0073 0.0070 0.0065 0.0058 0.0050 0.0041 0.0030 0.0020 0.0010 
0.90 0.0066 0.0065 0.0062 0.0058 0.0052 0.0044 0.0036 0.0027 0.0018 0.0009 
0.95 0.0058 0.0058 0.0055 0.0051 0.0046 0.0039 0.0032 0.0024 0.0016 0.0008 
1.00 0.0052 0.0051 0.0049 0.0045 0.0041 0.0035 0.0028 0.0021 0.0014 0.0007 
>1.00 0.057a 0.056a 0.054a 0.050a 0.045a 0.039a 0.031a 0.023a 0.015 0.007a 
az=e 2.4z 

The above functions have three important ly se 
characteristics: vaiomemee” 

I. In each case they consist of the boundary | 
function V,(1,2) plus other terms; hence the vil ly (Vi— Va) ky s-. S$ 
other terms may be regarded as a correction in as das ot- < on" , _— 
calculating the field at a general point. 

II. The first part of the correction consists of | —" Ss ” 
what might, by analogy with electric circuit ie "9 (19) 

( 


theory, be called a transient term. The term arises 
from the discontinuity in the boundary potential 
and decreases in a modified exponential manner 
as one leaves the region of the discontinuity. 

III. The second part of the correction, ap- 
plicable only for z>0 and #0, 1 consists of odd 
or even powers of z according as m is odd or even. 
The correction is necessary because V=2" does 
not satisfy Laplace’s equation except for n=0 
and n»=1. 

- As an example of the use of the above formulae 
in calculating potential fields consider an elec- 
trode system consisting of two coaxial cylinders, 
each of unit radius, separated a distance S along 
the axis (Fig. 1). The coordinate z=0 is chosen 
halfway between the ends of the cylinders; the 
left-hand cylinder is taken to be at a potential 
V,, the right-hand at a potential V». 

A first approximation to V(1,z) for this case 
can be taken as 


500 


in which it is assumed that the potential gradient 
in the gap between the cylinders is constant. 
Making use of the functions V,(1, 2) as defined 
in (14), V(1, z) can be expressed in the form 


V(1, 2) =Vo 
(Vi— Va) | 
+ —— 


r(net)-n(ve-3)| 


Thus the first approximation to the potential 
distribution becomes 


(20) 


Vir, 2) = Va 


(Vo— Va) | s . 
+————} V(r, +-) - v,(r, :--)}. (21) 
s | 2 2 


A better approximation is based on a Legendre 
polynomial expansion of the potential in the 
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gap. As is well known’ the Legendre polynomials 
P,(x) provide a series of functions which are 
orthogonal in the interval —1<x<1 and may be 
used to expand an arbitrary function in this 
interval. The P,(x) are polynomials of degree n. 
The expansion in terms of the polynomials has 
the valuable property of being the best least- 
square approximation to the function for a 
polynomial of degree m, where the series is 
terminated at P,,(x). 

If x in P,(x) is replaced by 2z/s the interval 
of orthogonality of the polynomials becomes 
—s/2<z<s/2; this interval is the one in which 
the expansion of V(1, z) is desired. The expansion 
takes the form 


~ 2 s s 
Vai, z) = > ¥ «,P,(-2), ~~ (22) 





n=0 AY 
where 
2n+1 /? 2 
a,=— f V1, 2)P.(-2)ads. (23) 
5 —s/2 Ss 


Since, for the particular case considered, V(1, z) 
is an odd function of z (excepting the constant 
term) and since P,,(x) is an odd or even function 
according as nm is odd or even, the integral in (23) 
vanishes for all even m (excepting n=0). As an 
analytical expression for V(1, z) is not known in 
the interval —s/2<z<s/2 the integral in (23) 
is evaluated numerically, the function V(1, z) 
having been determined in the interval by 
measurements made in an electrolytic plotting 
tank. The results of such measurements and a 
subsequent analysis for the coefficients a, are 


! 
Viz) Ve 





s 
¢ 
4 
Zz 


























Fic. 1. Boundary potential distribution of cylindrical 
electrode system. 


5L.R. Ford, Differential Equations (McGraw-Hill Book 
Company, New York, 1933), pp. 190-197. 
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TABLE VI. 
Separation S ao a1 a3 

R 0.500 0.38 0.054 
2R 0.500 0.35 0.049 
3R 0.500 0.35 0.052 
4R 0.500 0.36 0.052 
6R 0.500 0.37 0.057 
8R 0.500 0.37 0.044 








given in Table VI. For convenience the coeffi- 
cients are based on potentials Vz.=0 and V,=1. 
In terms of the cubic approximation the field 
at a general point can readily be found. Thus the 
potential V(1, z) can be approximated by 


Ss 
View s=-- 
2 
2 2 
Vat(Vi-— v.){ aoPo(-s) +a:P,( -2) 
b Ss 
V(1, 2) = (24) 
2 Ss Ss 
+a,P,(-s)], —-<2<-. 
Ss 2 2 
AY 
Vi, 2=-. 
2 





Since Po(x)=1, Pi(x)=x and P3(x) = }(5x°— 3x) 
the potential in the gap can be written 


2a, 20z* 32 
Vet(Vo- V.)|a0+—s-has —-=)| 
$s 


s8 S$ 


(2a; —3a3) 20a32 
s+ 





= Vat(Vi-—- V.){a0+ . (25) 


Ss s3 


Replacing z by [z+ (s/2) ]—s/2 and expanding in 
terms of [z+(s/2) ] one obtains 


(2a,+12a3) s 
Vat (Vi-— v.)| (ay~a1~as)4—— #+-) 


5 


30az s\? 20a; s\3 
02) -(02)| 
s? 2 s3 2 


so that in terms of the functions V,(1,2) the 
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gap potential can be expressed 
s 

Vat (Ve— v)| (do — 4, — 43) vo( z+ ) 
2 


(2a; + 12a 3) Ss 30a 3 Ss 
+ v,(1, 2+ )- v(1,2+-) 
s 2 s* 2 
20a; 


ay 
V(1, 2+ ) . (26) 
s3 2 


Similarly, by writing [z—(s/2)]+s/2 for z the 
corresponding expression becomes 


+ 


s 
Ve" vo( 1, 2) +(Vs- vo (ao +a,+as3) 


Ss (2a,+12a3) s 
tia Od 02) 
2 Ss 2 
30a 3 | 9 AY 20a 3 J 3 Ss 
+ (1,5-=)+- -— (1.2 )}. (27) 
s? 2 s3 2 


The expression (26) represents the potential 
V(1, 2) for z<s/2 since for z< —s/2 it reduces 
to V,. For z<s/2 (27) is zero while for z>5s/2 it is 
identically equal to (26). Thus subtracting (27) 
from (26) and adding V,-Vol1,2—(s/2)] the 
desired form of V(1, z) is obtained. The potential 
at a general point is thus approximated by the 
formula 


V(r, 6) = Vert(Vo= Va){ (aaa) 
s s 
| Vol ret) + ro( :--)| 
2 2 
(2a,;+12a3) Ss Ss 
+ AL vi(1. 242) V(r, 2-2) | 
Ss 2 2 
30asf AY s\) 
amen v(r,24-)+02(r,5--) 
s? 4 2 2 


20a3f s S\ 1 
oo — v:(r.2+-)-Va(r.2--) ' (28) 
et 2 2/ 3 


With the aid of Tables II-V the potential 
field for the electrode system under discussion 
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Fic. 2. Equipotential linés of cylindrical electrode system 
(V.=0,V,=8, s=0.2R). 


can now be calculated. The equipotential lines 
for the case s=0.2R are plotted in Fig. 2. These 
were obtained from a calculation using the linear 
approximation, taking @9=0.500, a, =0.350, and 
a3;=0.000. For comparison, data obtained in an 
electrolytic plotting tank are shown as circled 
points. The discrepancies can be attributed to 
inaccuracies of the measurements. The error 
arising from the linear approximation is greatest 
in the vicinity of the edges of the electrodes. 
If greater accuracy is required the higher order 
approximation may be used. 

The extension of the method to more general 
cases of axially symmetric electrode systems 
having any electrode configuration is obvious. 
The only requirement is that the potential along 
a cylinder of constant radius—within which 
Laplace’s equation holds—be known with the 
required accuracy. 

The essential principle of the method outlined 
in this paper for the calculation of fields is that 
of superposition. Complex boundary potential 
functions are expressed over short intervals in 
terms of functions which can be evaluated and 
the results superposed to obtain the complete 
solution. Tabulated values have been included 
to facilitate the use of the method. 
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The Operational Theory of Longitudinal Impact 


Louis A. PIPES 
Graduate School of Engineering, Harvard University, Cambridge, Massachusetts 


(Received June 2, 1942) 


The differential equations governing the longitudinal displacement and velocity of a uniform 
rod are identical to those governing the propagation of charge and current along a uniform 
electrical transmission line. By making use of this analogy hetween the rod and the transmission 
line, it is possible to translate some well-known results of electrical circuit theory to mechanical 
terms. In this manner, and by the use of the operational calculus, interesting and useful results 
in the theory of longitudinal impact are obtained. The treatment is vastly simpler and more 
general than the classical theory of Saint-Venant and Boussinesq and some of the results are 


thouglit to be new. 





I. INTRODUCTION 


HE problem considered is that of determin- 
ing the subsequent behavior of a long, 
uniform thin rod that is struck at one end by a 
massive body moving in the direction of the 
length of the rod. Different states of mobility 
of the other end of the rod are assumed. The case 
where the rod is embedded in a viscous medium 
is also studied. 

The classical solution of this problem was ob- 
tained by J. Boussinesq'? and Saint-Venant® 
using the general solution of the wave equation 
that governs the propagation of the displacement 
of the rod. This method of solution is laborious 
and difficult to extend to general cases of interest. 
The method presented in this paper is much 
simpler and permits more general cases to be 
studied than is possible by the method of 
Boussinesq and Saint-Venant. 


Il. THE EQUATIONS OF MOTION- 
ELECTRICAL ANALOGY 


Consider a uniform rod or bar as shown in 
Fig. 1. 





L J 





om 
—_— > 


Fic. 1. Uniform rod or bar. 


1J. Boussinesq, Applications des Potentie's a l'étude de 
Véquilibre et du. Mouvement des Solides Elastiques (Paris 
1885), p. 508. 

2A. E. H. Love, A Treatise on the Mathematical Theory of 
Elasticity (Cambridge, 1934), p. 431. 

’ Barre de Saint-Venant, Theorie de l'elasticité des corps 
solides (Paris, 1883), p. 307. 
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Let us introduce the following notation: 


E,=the Young’s modulus of the material of 


the bar. 
p=the density of the material of the bar. 
t=the time. 


x=the longitudinal coordinate measured 
from one end of the bar. 

u(x, t)=the longitudinal displacement of any 
cross section of the bar during vibration 
from its position of equilibrium. 

v(x, ) =the velocity of any cross section of the 
bar during vibration. 

F(x, t)=the longitudinal compressive force act- 
ing across any section for the bar. 

A =the cross-sectional area of the bar. 


It is easy to show‘ that the displacement, 


‘ u(x, t) satisfies the wave equation, 


0?u / dt? = Co?d?u/dx? (1) 
co*=E,/p (2) 


where 


and x is measured as in Fig. 1. The compressive 
force F(x, t) satisfies the equation, 


F=—AE,du/ox. (3) 


By differentiating (1) with respect to x and sub- 
stituting (3) into the resulting equation, we 
obtain: 

82°F /dt? =c2d2F/dx?. (4) 


That is, the force F(x, t) satisfies the same 
wave equation that is satisfied by the displace- 
ment, u(x, ft). The velocity, v(x, t) of any point 


4S. Timoshenko, Vibration Problems in Engineering 
(D. Van Nostrand, 1937). 
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of the rod is defined by, 
v(x, t) =du/dt. (5) 


Differentiating (1) with respect to ¢ and sub- 
stituting from (5), we obtain, 
070 / Ol? = Co?0*v/dx?. 


(6) 


Differentiating, (3) with respect to t, we have 


oF 0 0 
—=—(AE,)—(du/dx) = —(AE,)—(0u/0l) 
ot ot Ox 
= —(AE,)(dv/dx). (7) 
Equation (1) may be written in the form, 
pd*u/dt? = E,d*u/dx?. (8) 
Substituting (3) and (5) into (8), we obtain, 
ov 1 OF 
p—=-——. (9) 
at A 0x 


Equations (7) and (9) may be written in the 
form, 
1 oF ov 
—=——, (10) 
Ox 


pAdv/dt= —dF/dx. (11) 


Let us now consider a dissipationless electrical 
transmission line as shown in Fig. 2. 








x — 





Fic. 2. Electrical transmission line. 


Let L=the inductance per unit length of the 
transmission line. 
C=the capacitance per unit length of the 
transmission line. 
e(x, t)=the potential difference between the line 
conductors. 
i(x, t) =the line current. 
x =the space coordinate measured from one 
end of the line. 
t=the time. 


It may be shown’ that the equations governing 
the propagation of potential and current along 





5 E. A. Guillemin, Communication Networks (John Wiley, 
1935), Vol. 2, Chap. 1. 
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the line are, 
Cde/dt = —0d1/dx, (12) 


Lai/at= —de/ax. (13) 

The Eqs. (10) and (11) are of exactly the same 
form as (12) and (13). We thus see that we may 
write the following relations, 


F(x, t)= e(x, t), (14) 
v(x, t)=14(x, 2), (15) 
1/AE,=C, (16) 
(pA)=m=L. (17) 


The symbol + expresses the fact that the two 
quantities are analogous and m is the mass per 
unit length of the rod. By the above analogy the 
mechanical problem may be translated into an 
electrical problem. Many well-known results in 
the theory of electrical networks may be used. 
To do this, it is only necessary to replace masses 
by inductances, springs by suitable condensers, 
and frictional forces of the coulomb type by 
electrical resistances. 

The problem of determining the behavior of a 
rod of length s and fixed at one end when struck 
at the other end by a mass mp moving with a 
velocity vo as shown in Fig. 3, is analogous to the 





ANANANS 





SE] 


Fic. 3. Fixed bar. 











\ 
te — 


x=0 x=8 





Fic. 4. Electrical equivalent of fixed bar. 


problem of determining the current and poten- 
tial distribution of an open-ended transmission 
line when an inductance Lo carrying an initial 
current 7%» is suddenly connected to the other end 
of the line. In order to facilitate the translation 
of electrical to mechanical quantities, the fol- 
lowing table of analogs, is useful. 
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TABLE I. Table of analogs. 


Rod Transmission line 





displacement, u(x, ¢) 

velocity, v(x, t) / 

compressive force, f(x, t) 

lumped mass, mo 

distributed mass/length 
m= (pA) 

compliance/length, 1/AE, 

velocity of wave propaga- 
tion, ¢o= (E,/p)* 

initial velocity, vo 

length or rod, s 


charge distribution, q(x, t) 
current, 7(x, ¢) 

potential distribution, e(x, t) 
lumped inductance, Lo 
inductance/length, L 


capacitance/length, C 

velocity of wave propaga- 
tion, %=1/(LC)! 

initial current, Zo 

length of transmission line, s 

ratio of striking mass to ratio of inductance Lo to 
mass of rod, m= mo/ms line inductance n= Lo/sL 

damping coefficient per unit _resistance/length, R 
length, Rn 

dashpot damping coeffi- 
cient, Ra 

spring stiffness, ko 

mechanical surge imped- 
ance of rod, com 


resistance, Ro 


elastance, So 
surge impedance, 
Zo=(L/C)* 


Ill. THE OPERATIONAL EQUATIONS OF 
THE TRANSMISSION LINE 

By the above analogs, the study of longi- 
tudinal impact of rods or bars is reduced to the 
study of transients along a transmission line 
which are initiated by certain initial conditions. 
The simplest mathematical technique to use in 
solving problems of this character is the opera- 
tional calculus. 

As is well known, the modern theory of the 
operational calculus is based upon the Fourier- 
Mellin integral equations.® This set of equations 
states that if, 


g(p)=pf e*h(ddr, (18) 
then, : 
1 c+j@w 
h(t) -— |{ g(p)e”‘dp/p, where j= y/ —1 (19) 
2rj c—jo 


and ¢ is a real constant of such magnitude that 
the poles of the integrand lie entirely to the left 
of the path of integration in the complex p plane. 
For brevity, let us use the notation, 


g(p) =Lh(t), (20) 

h(t)=L~'g(p), (21) 

to express, symbolically, the relations (18) and 
(19). 

Let us now consider the transmission line of 


Fig. 5. 


6° L. A. Pipes, ‘The operational calculus,” J. App. Phys. 
10, 172, 258, 301 (1939). 
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Z, and Z, represent the impedances of linear 
networks connected to the sending and receiving 
ends of the line. At the sending end, an electro- 
motive force e,(t) which may be a general func- 


4 (x, 4) —» 


ee ol | [z.] 


A(x 4) 











tL Lat 


xX=o ate 








Fic. 5. Transmission line with general impedances. 


tion of time is impressed. If the line is dissipa- 
tionless and has an inductance L and capacitance 
C per unit length, then the potential e(x, ¢) and 
the current i(x, ¢) satisfy the Eqs. (12) and (13) 
along the line. The equations 


e,=Z,1(0, t)+e(0, 2) (22) 
e(s, t)=Z,1(s, t) (23) 


and 


must be satisfied at the sending and receiving 
ends, respectively. If we now introduce the 
transforms, 


E,(p) =Le,(t), (24) 
E(x, p) =Le(x, t), (25) 
I(x, p) =Li(x, t), (26) 


and make use of Eqs. (12), (13), (22), and (23) 
we obtain,® 
; ZoE,[e**-” +7r,e-*¢-” | 
E(x, p)= —, (27) 
(Z.+Zo)(e®—r re) 


E,[e*“s-» —r,e-os-2) | 








I(x, p) = ; (28) 
(Z.+Zo) (e* —r,r,e-**) 
where 
a=p(LC)}, (29) 
Zo=(L/C)}, (30) 
°,= (Z,—Zo) /(Zs+Zo), (31) 
’,= (Z,—Zo)/(Z-+Zo). (32) 


These are the well-known equations for the 
transforms of the potential E(x, p) and current 
I(x, p), respectively, of a dissipationless line. 
The quantities r, and 7, are the operational 
reflection coefficients and Zo is the characteristic 
or surge impedance of the line. 

If the character of the impressed electromotive 
force e,(t) is known, then it is possible to com- 
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pute its transform E,(p) and the potential and 
current distribution may be determined by 
finding the inverse transforms of E(x, p) and 
I(x, p) as given by (27) and (28). In the theory of 
impact, we are in general interested in the so- 
called wave solution of the problem. That is, 
the impact of the mass mp on the rod produces a 
wave of force and deflection that travels along 
the rod with a certain velocity and is reflected 
at the end of the rod in a certain manner. To ob- 
tain the character of these waves, let us divide 
the numerator and denominator of (27) and (28) 
by e*. We then obtain, 


ZoE,[e-* +r ,e-8—-» 7 


(Zo+Z,)(A— rare) 





E(x, p)= (33) 


E,[e-™ —_ 7g e(2e—2)] 


I(x, p)= (34) 


7 (Zo+Z,)(1 ie rere) 


Ihe second factor in the denominator of (33) 
and (34) may be expanded in the form, 


w® 
[1—r,r,e-2* =D (r.r,) "e-2*. (35) 
n=0 


Substituting this into (33) and (34), we obtain, 


LoE, x 
E(x, p)= me DS (rr) eat) 
; (Z.+Zo) n=0 
ZoE,r, w 


eae Zz (r fr) mg —e(2e+2en—z) 


(36) 


.. » 
cinema > (r.%,) neg—a(z+2sn) 


(Z,+Zo) »=0 


I(x, p)= 


E,r, «@ 
2 (r.7,) Ne—a(2s+2sn—z) 


(Z,+Zo) n=0 





(37) 


The expressions (36) and (37) have simple 
physical interpretations. First let us consider 
the case where the terminal impedance of the 
system Z, satisfies the equation 


Z,=Zo=(L/C)'. (38) 
In this case, we have from (32), 
r,-=0. (39) 





In such a case (36) and (37) reduce to the simple 
forms, 


E(x, p) =ZoE,e-**/(Z,+Zo), (40) 
I(x, p) = E,e~**/(Z,+Zo). (41) 


These expressions are the transforms of waves 
of potential and current progressing along the 
line with a velocity of, 


Vo=1/(LC)}. (42) 


These waves are not reflected from the re- 
ceiving end as a consequence of the relation 
r-=0. If r,+0, the series of terms in (36) and 
(37) are the transforms of waves produced by 
successive reflections from the ends of the line. 


IV. ROD FIXED AT ONE END AND STRUCK 
LONGITUDINALLY AT THE OTHER 


Let us consider the case of Fig. 3. In this case, 
the rod is struck by a mass mo moving with 
velocity vo and the other end of the rod is fixed. 
To solve the problem of determining the sub- 
sequent state of affairs, we consider the electrical 
analog of the rod as shown in Fig. 4. 

In order to use the general results (36) and 
(37) for the case of the inductance Ly carrying a 
current %) and suddenly switched on to the open- 
circuited line, we place, 


E,=Lopio, (43) 
Z,=Lop, (44) 
Z,=@, (45) 
Hence we have, 

rr=(Z,—Zo)/(Z-+Zo) = 1, (46) 
rs=(Lop—Zo)/(Lop+Zo), (47) 

E,/(Zs+Zo) =Loplo/(Lop+Zo) 
=pIo/(b+Zo/Lo). (48) 


Then by (36), we have: 


Zopto 
E(x, p) ee 
P+Zo/Lo 





[e-** + 7 e-4(28—2) 


+r re-use) +r. ees») + ewes ] (49) 


for the transform of the potential. 
To interpret this result, let us consider the 
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first term of (49). 


Zopive—** Zopto exp ( —xp/ Vo) 





: piginnmuatemchnems,. iit 
(p+Zo/Lo) (p+) 
where, 
b=Zy ‘Lo, (51) 
Vo=1/(LC)}. (52) 


From the theory of Laplacian transforms,® we 
have the following inverse transform, 


__Zobio exp (—bx/ V0) 
(p+) 
0 ft <x/ Vo 





4 


lieZo exp (—Db(t—x/Vo)) t>x/Vo. (53) 


The first term of (49) thus represents the 
inverse transform of a traveling wave of electro- 
motive force which progresses along the line with 
a velocity of Vo. It is convenient to introduce the 
quantity, 

Ly/Ls=n. (54) 
n thus represents the ratio of the lumped in- 


ductance Ly to the entire inductance of the line. 
With this notation, we may write, 


ipZo exp [ —b(t—x/ Vo) ] 
=i9Z 9 exp (x—tVo)/ns]. (55) 
To interpret the mechanical significance of the 
solution, it is necessary only to examine Table I. 
From this, we have, 
ipZy exp [—(t—x/Vo)b ] 


(0 t<x/co 
= { _ (56) 

\vocom exp [(x—cot)/ns] t>x/c. 
That is, the effect of the impact is to create a pres- 
sure wave which progresses along the rod to the 
fixed end with a velocity co. Let us now examine 
the inverse transform of the second term of (49). 
Since r, is equal to 1, this term is, 


Zopioe*°*—”) 
r —1 





4 


(p+5) 
(0 t<(2s—x)/Vo 
= (ipZo exp [ —b[t—(2s—x)/Vo]] (57) 
| t>(2s—x)/Vo. 
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By Table I, we have: 
L-%49Z exp [—bd[t — (2s —x)/ Vo) 
= vocom exp [[(2s—x)—cot ]/ns]. (58) 


We thus see that the second term of (49) is the 
transform of a wave of compression traveling 
back from the fixed end of the rod with a velocity 
co. In the same manner, it may be demonstrated 
that the inverse transform of the third term of 
(49) is given by, 


Zopior e742) 
(p+b) 
0 t<(2s+x)/Vo 
ioZo exp [—b[t—(2s+x)/Vo]] 





a= , _ 2ioZ ob[ t mat (2s+x)/ Vo] (59) 
| Xexp | —bLt—(2s+x)/Vo]] 
{ t>(2s+x)/Vo. 


By Table I, we have, 
ioZo exp [ —b[t—(2s+x)/Vo ] 
— 2ipZoblt —(2s+x)/ Vo] 
Xexp [—bLt—(2s+x)/Vo}] 
= vocom exp [[(2s+x) —cot ]/ns | 
X {1+2[(2s+x) —cot]/ns}. (60) 


This represents a wave of compression travel- 
ing from the end of the rod that was struck 
toward the fixed end with a velocity co. If the 
striking body becomes attached to the rod, we 
may calculate the subsequent waves of pressure 
by determining the inverse transforms of the 
various terms of (49). However, if the bodies 
remain detached, the solution continues to hold 
so long as there is positive pressure between the 
rod and the striking body. 

To examine when the pressure ceases to be 
positive, we place x=0 in the wave solutions. 
If we let fo be the compressive force at x=0, 
we have, 


fo=vocom exp (—Cot/ns) 
for 0<t<2s/co. (61) 


This expression is positive. However, in the 
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range, 2s/cg<t<4s/co, the compressive force is 
given by, 


fo=vocom exp [ —cot/(ns) ] 


Col _ 2s 
x| 1+2e(1—" — )| (62) 
ns 
This expression vanishes, when, 
2cot/ns =4/n+2+e-?/". (63) 


It is possible for this equation to have a root 
in the interval 2s/cy9<t<4s/co if, 


2+e77/"<4/n. (64) 
It may be shown that the equation, 
2+e-2/"=4/n (65) 


has a root lying between n=1 and n=1.73. 
Hence if »<1.73 the impact ceases in the above 
interval and the time at which it ceases is given 
by, 


t=s(2+n-+ne-?!"/2)/co. (66) 


To interpret the wave solution for the velocity, 
it is necessary only to compute the inverse 
transforms of I(x, p) when the values given by 
(43), (44) and (45) are substituted into (37). 


V. ROD FREE AT ONE END AND STRUCK 
LONGITUDINALLY AT THE OTHER 


Let us now consider the case of Fig. 6. In this 


*~EC 7 


Fic. 6. Free bar. 








case, the bar is free to move and the electrical 
analog of the system is an inductance con- 
nected to a short-circuited line at t=0, as shown 


in Fig. 7. 


Fic. 7. Short-circuited line. 











In order to treat this case, it is necessary only 
to realize that we now have, 


Z,=0. 


(67) 








and hence, 


r,-=—1. (68) 
The relations (43), and (44) still hold. From 
(49), we have, on placing x=0 and 7,= —1. 
E@, =e 1 
1( _?~p ine emnbsindiiea —e 2as__y ¢ 2as 
b) 
et... J, (69) 


This is the transform of the potential at x=0. 
In view of the transforms (53), (57), and (59) 
it is easily seen that the potential, and hence the 
mechanical compressive force becomes negative 
at x=0 at the instant, 


t=2s Co. (70) 


This means that the impact produces a wave of 
compression which is reflected from the free end 
as a wave of extension that causes the impact 
to cease. It is interesting to compute the ve- 
locity of the rod after impact. 

To do this, we substitute (43), (44), (47), and 
(68) into the transform of the current given by 
(37) and obtain, 





pro 
I(x, pb) Bical a teiealeias [e-*7 + e222) 
(p+) ; 
—re ost...) (71) 
From (53), we have 
pice~** (0 t<x/ Vo 
L-'— ={ (72) 


(p+b) | ige(=—V ot) / ns t>x/Vo. 


Hence the mechanical analog of the first term 
of (71) is the transform of a wave of velocity 
given by (72). From (57), we have, 

; 0 t<(2s—x) Vo 
pie 2'2#—) ] 


-———— =igexp[[(2s—x)—Vot]/ns] (73) 


| 
| 
( t>(2s—x) Vo. 


The mechanical analog of the right member 
of (73) represents a wave of velocity reflected 
from the free end and traveling toward the struck 
end with a velocity co. At t9=2s/co the velocity 
of the rod is given by, 


v(x, to) = Vo[ e~2/ "e2! net e—z/ns'] (74) 


or if we let 
x=sS—yY, 


(75) 
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this reduces to 


2ve—"/" cosh (y/ns). (76) 


It is thus obvious that the rod rebounds vi- 
brating. 


VI. ROD TERMINATED BY A DASH POT 


A very interesting case arises when the rod is 
struck at one end and the other end is connected 
to a dash pot of a certain strength, as shown in 
Fig. 8. 





+AC _ 





Fic. 8. Rod terminated by a dash pot. 
The electrical analog of this is the transmis- 


sion line terminated by a resistance as shown 
in Fig. 9. 


] 


Fic. 9. Line terminated by a resistance. 
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In this case, we have, 


Z,=Ro (77) 
r,=(Ro—Zo) ‘(Ro+Zo). (78) 


and 


It is thus seen that the reflection coefficient 7, is 
a numerical factor less than one. Let us consider 
the case where the electrical transmission line is 
terminated in a resistance whose value is equal 
to the characteristic impedance of the line. 
That is, the case where 


Ro=Zo. (79) 


In the mechanical analog, this means that 
one end of the rod encounters a resistance pro- 
portional to its velocity of magnitude, 


Rn=com. (80) 

In this case, we have, 
7r,=0, (81) 
and the solution is vastly simplified because we 
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do not have any waves reflected from the dash- 
pot end. In this case, the transforms of the 
potential and current are given by, 


E(x, p) = Zoptoe~**/(b+5), (82) 
I(x, p) = pice**/(p+b). (83) 


The transform of the charge distribution is 
given by, 


Q(x, p) =I (x, p)/p=toe**/(p+b). (84) 


Interpreting the inverse transforms of these 
quantities in terms of their mechanical analogues, 
we obtain the following expressions for the 
compressive force, velocity, and displacement of 
the rod, 


{° t<x/Co 

f(x, )= an (85) 

\vocom exp [(x—cot)/ns] t>x/co, 
t<x/co 

v(x, t) = (86) 
\voexp [(x—cot)/ns] t>x/co, 
'y t<x/Co 

u(x, t)= [1—exp (x—cot)/ns luosn/co (87) 
| t>x/Co. 


From (87) we have the result that the dis- 
placement of the struck end is given by . 


u(0, t)=[1—exp (—cot/ns) Juosn/co. (88) 


It is thus evident that as a consequence of the 
impact, the struck end moves a distance snvo/Co. 


VII. INFINITE ROD UNDER ACTION 
OF A SPRING 


To illustrate how the general equations may 
be used to analyze special cases, let us consider 
the case of Fig. 10. In this case, a spring of 





}-4-- 


3 








am, 


Fic. 10. Spring and infinite rod. 
spring constant ko is deflected from its position of 


equilibrium by a distance wu» so that it has an 
initial potential energy of kouo?/2. The spring 
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is released at t=0 and its potential energy is 
transferred to an infinite rod. 

The electrical analog of this system is a 
condenser discharging into an infinite trans- 
mission line as shown in Fig. 11. 








Fic. 11. Condenser discharging into infinite line. 


In the electrical analog, the electric energy 
of the condenser is SoQo?/2, where So is the 
elastance of the condenser and Q2, is its initial 
charge. In order to determine the waves of 
potential and current we place, 

r,=0, 


E, = SoQo, 
Z; = So p, 


(88) 
(89) 
(90) 
in the general Eqs. (36) and (37). We thus obtain, 
(Zo+So/p), 
I(x, p) = SoQoe~**/(Zot+So/p), 
Q(x, p) =I(x, p)/p. 


E(x, bp) = SoQoZ ve?” (91) 


(92) 
(93) 


Computing the inverse transforms of these 
expressions, we obtain, 


(0 t<x/Vo 
e(x, t)=4 


94 
| SoQoe Soe (2-t¥ 0) (94) 


t>x/Vo, 
(0 t<x/Vo 

i(x, t) = ( SoQo ie : 
|\—— exp [SoC(x—tVo)] t>x/Vo, (95) 
L Zo 

‘0 t<x/Vo 


q(x, t) =: (96) 
Qo[1 —exp [ SoC(x —tVo) | t>x, Vo. 


These expressions may be easily translated 


into mechanical language by the table of 
analogs. 


VIII. LONG ROD EMBEDDED IN A 
VISCOUS MEDIUM 


A problem of some interest, is that of de- 
termining the character of the waves propagated 
along a rod that is embedded in a medium that 
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impedes its motion with a force proportional 
to its velocity. To simplify the problem, let us 
consider the case of an infinite rod or one that 
is terminated by a dash pot of such a strength 
that no waves are reflected. We shall consider 
the case in which a force fo is suddenly applied 
to the rod as shown in Fig. 12. 


@ dbf LALLA 
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Fic. 12. Rod embedded in a viscous medium. 


The electrical equivalent of this system is 
given in Fig. 13. 





& 





Fic. 13. Electrical transmission line with resistance. 


In this case, the electrical transmission line 
has a resistance R per unit length, and hence the 
operational propagation constant, a, is given by 


a=[(Lp+R)(Cp) } 


instead of by (29). 
To determine the waves of potential and cur- 
rent, we place, 


(97) 


r,=0, (98) 


E,=Eo, (99) 


Z.=0, (100) 


in the general Eqs. (36), and (37). We then have, 


E(x, p) = Eye~*?, (101) 


I(x,-p) = Eve-**/Zo. (102) 


In this case the computation of the inverse 
transforms of (101) and (102) is a matter of 
some complexity, because of the complexity of 
the operational propagation constant a and the 
operational characteristic impedance Zp») which 
is now given by, 


Zo=[(Lp+R)/Cp}}. 


The required inverse transforms may be 
computed by the use of the basic contour in- 
tegral of Eq. (19). To facilitate the writing, 


(103) 
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let us place, 
a=R/2L. (104) 
We then have, 


a=(LC)*|p(p+R/L) }}=[p(p+2a)/Vo]}'. (105) 
The inverse transform of (101) is then given by,’ 


0 t<x/Vo 


Bal 1" 


Ee~*? = ¢ (106) 
ax ct e-*D,[a(u?—x? oe] 


Vo 2z/Vo (u? —x?/ Vo?)} 





t>x/ Vo, 


where J,(y) is the modified Bessel function of the 
first order and argument y. The transform (106) 
shows that a wave of potential or a wave of 
compressive force in the mechanical case ar- 
rives at a point x from the transmitting end at 
a time fj=x/Vo. At that instant, the potential 
or compressive force rises suddenly to the 
value Eo exp (—x/Vo). The compressive force 
wave thus travels along the rod with a vertical 
front whose height steadily decreases with an 
increase in x. As time progresses, the potential 
assumes the value Ey throughout. 

To compute the inverse transform of the 
current, we write, 


I(x, p) => Eve—**/Zo 
Eo(C/L) ‘p exp L—x[ (p+ 2a) p |*/ Vo] 


- [(p+2a)p]} 


7N. W. McLachlan, Complex Variable and Operational 
Calculus (Cambridge, 1939), pp. 216-218. 





(107) 
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The inverse transform of this expression is 
given by,’ 


0 t<x/Vo 
L— I(x, p) = / Eo(C/L)'e-#*I (108) 
| X[a(t?—x?/Vo?)#] t>x/Vo. 


From (108) it is seen that in the analogous 
mechanical case, a wave of velocity arrives at a 
distance x from the transmitting end at a time 
to =x/Vo. At this time the current rises suddenly 
to the value E[ (C/L) ]* exp (—ax/Vo). As time 
progresses, the current finally assumes the value 
0. To obtain the mechanical displacement of the 
rod, it is necessary to compute the analogous 
quantity, 


q(x, =f i(x, t)dt. (109) 


This expression gives the wave of displacement 
that progresses along the rod when interpreted 
in terms of the various mechanical analogs. 


IX. CONCLUSION 


It is apparent from the above examples that 
treating the problems of longitudinal impact of 
rods in terms of their electrical analogs vastly 
simplifies the analysis. More complicated termi- 
nal conditions than those treated in this dis- 
cussion may be analyzed by introducing proper 
terminal impedances in the basic Eqs. (36) and 
(37). The classical theory of Saint-Venant and 
Boussinesq is difficult to extend to general cases. 
The analysis presented in this paper emphasizes 
the utility and simplicity of operational methods 
in problems dealing with impulsive forces. 





Energy and Energy Flow in the Electromagnetic Field 


J. SLEPIAN 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 
(Received April 17, 1942) 


A new method for deriving the usual Poynting vector is given. This method also yields other 
equally valid energy flow vectors with the same or other equally valid postulated electro- 
magnetic energy densities. Examples of such alternative Poynting vectors with their associated 
energy densities are given. A definition is given which distinguishes between the conduction and 
displacement currents in matter. It is shown that by addition of a term, the postulated energy 
flow widely used by electrical power engineers and based on wattmeter readings becomes a 


valid alternative Poynting vector. 


] Poynting’s theory of the energy streaming 
¢ in an electromagnetic field is usually based 
on the interpretation of the various terms in the 
following equation, derived directly from Max- 
well’s equations: 


Cc 
—div [EXH] 


T 


c c 
= ——/]-curl E+—E-curl H 


4r 4r 
c 1 0B 
‘ u-(-—)+e-1 
4n c él 
1 0B 
=—f]. +E-I;. (1) 
4n ot 


In this equation, J; is the total current density, 
i.e., the sum of the conduction, convection, and 
displacement current densities, and perhaps best 
defined as (c/42) curl H. The other vectors have 
the significance usually attached to them in 
Maxwell's theory. 

By examining the meaning of the terms on the 
right-hand side of this equation in various special 
cases, it may be made plausible that these terms 
represent the rate of increase of an electromag- 
netic energy density to be assigned to the point 
in the empty or material-occupied space being 
considered, plus the rates of increase of other 
recognizably localized energy densities such as 
heat, chemical energy, etc., which are to be 
regarded as due to the electromagnetic field, plus 
the rate of doing work of certain ponderomotive 
forces acting on moving matter per unit volume 
at the point in question. These ponderomotive 
forces are regarded as arising from, or being due 
to, the action of the electromagnetic field on 
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charged, electrically polarized, or magnetized 
matter. With this interpretation of the right- 
hand side of this equation, an appropriate mean- 
ing may be given to the Poynting vector, 
(c/4r)_E XH], which appears on the left-hand 
side of the equation by postulating that it 
represents a flow of ‘electromagnetic’ energy in 
space or matter. 

It is well known that this Poynting vector is 
not uniquely defined by Eq. (1), but any other 
vector with an equal divergence can equally well, 
according to Eq. (1), claim to represent a flow 
of electromagnetic energy. The method of de- 
riving energy flow vectors given in this paper 
exhibits these alternative Poynting vectors with 
a kind of physical interpretation. 

It is probably fairly evident, also, that the 
success of the Poynting vector in Eq. (1) depends 
on the choice of the electromagnetic energy 
density distributed through space. It is also 
clear that this choice of energy density is to a 
considerable extent arbitrary, since only its total 
quantity, integrated over all space may be said 
to be observable, so that any other assumed 
density, so long as it gives the same total when 
integrated over all space, may be regarded as 
equally valid. Of course, such other assumed 
electromagnetic energy density will have corre- 
sponding to it a valid energy-flow vector which 
will differ from the Poynting vector even in its 
divergence. This paper shows how such other 
energy flow vectors may be constructed. 

2. The kinds of special cases which are con- 
sidered as establishing the interpretation of the 
terms in (1) in general can be reduced to cases 
where the changes in the electromagnetic field 
are produced by the motion of electric charges 
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alone. By this we mean that energy of a mechani- 


cal form changes in magnitude at a rate given 


by the summation of pE-v over all the charged 
matter, where p is the charge density, E the 
electric field intensity, and v the velocity of the 
charged matter, and that no changes in any 
other recognizable form of energy can be dis- 
cerned, other than the postulated electromag- 
netic energy which is being calculated and the 
appearance of heat which must be removed if 
the temperature at any point is to remain un- 
changed. The rate of change of the postulated 
electromagnetic energy is set equal to the rate of 
increase in mechanical energy due to the motion 
of the charges plus the rate of development of 
heat. The formulae for the electromagnetic 
energy change calculated for these special cases 
are then assumed to be universally valid. 

Further independent assumptions are then 
made as to the distribution of this energy in 
space. The formulae derived in the special cases 
are, or may be converted into, integrals over 
space of functions of the field vectors, i.e., E, D, 
H, and B. These functions of the field vectors are 
then postulated to be the electromagnetic energy 
density in the special cases considered, and also, 
so far as they may be given meaning, in general 
cases. Where the formulae derived for the special 
cases are not applicable, the Poynting vector is 
used to define an electromagnetic energy density 
by postulating that an energy density function 
T exists defined by 


6T 6Q _ € 

—-+—=-div —[EXH], (2) 

ét = bt 4a 
where 6/6tis the convective derivative for moving 
matter, and 6Q/ét is the rate of appearance of 
heat per unit volume of matter. Thus, for ex- 
ample, Eq. (2) for magnetic matter, in general, 
makes the energy density, J, be given by 


5Tn 60 1 6B 
—+—=—H-—. (3) 
bt ét 4 ét 
Similarly for dielectric matter, the electric energy 
density is defined by 
6T. 60 


1 6D 
+—=E-I,=—E-—+E-I,, (4) 
) 4nr = 6t 





where J, is the conduction current density. 


VOLUME 13, AUGUST, 1942 


3. In this paper we consider not special cases, 
but proceed at once to the most general electro- 
magnetic field. However, for simplicity, we shall 
consider only fields with matter at rest, although 
the method is readily extended to matter in 
motion. We take an arbitrary volume of this 
electromagnetic field, enclose it in a geometric 
closed surface, and define an energy flow through 
this surface in the following way. 

We suddenly create a narrow empty shell 
following the surface. Within this shell we place 
means for observing the defined energy flow. 
Thus the definition we are developing for energy 
flow uses a procedure similar to that frequently 
used for defining D and B within matter. 

The requirement for ‘‘suddenness”’ in the crea- 
tion of the empty shell is for the purpose of 
distinguishing between the displacement and con- 
duction current. In such a “‘suddenly”’ created 
shell the normal component of electric force will 
equal the normal component of D. The rate of 
change of the normal component of electric 
force will be 4x times J;, the total current 
density. 

In the empty shell we place an inner shell of 
magnetic material of infinite permeability. This 
magnetic shell we expect to cut off all interplay 
or interchange of energy between the electro- 
magnetic field inside the shell and the field 
outside. 

To compensate for the disturbing effect of the 
magnetic shell on the inner field, we place 
charges between the shell and the field in proper 
amounts, and give them such motions that the 
field inside is the same and for some time con- 
tinues to be the same as if the shell space had 
not been created. If we like, we may also place 
moving charges outside the magnetic shell to 
keep the external field unchanged. 

The rate at which work is done by the moving 
charges is then defined as the energy flow out 
through the surface from the inner field, with due 
consideration of possible changes in surface 
energy density at the shell itself. 

The creation of the empty shell, because of its 
infinitesimal thickness, may perhaps be regarded 
as involving a negligible amount of energy in its 
volume. However, where the shell passes through 
matter, discontinuities are produced in the field 
vectors, and the possibility must be considered 
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that a surface distribution of energy must be 
assigned to these discontinuities. 

Likewise, the bringing in of the infinitely 
permeable material to make the inner magnetic 
shell, and the bringing in of the compensating 
charges, presumably from infinitely remote re- 
gions, certainly required the expenditure of work, 
and therefore the storage of a certain amount of 
energy which may perhaps need to be ascribed to 
the shell. However, this energy of creation of the 
shell with compensating charges does not come 
into the picture because only the changes which 
occur after the creation of the shell are considered 
in the calculations. 

But, in moving the compensating charges 
about, local microscopic electric and magnetic 
fields about them will change, and it is possible 
that a changing surface distribution of energy 
will need to be assigned to these moving com- 
' pensating charges. 

It turns out that equally self-consistent but 
different Poynting vectors are obtained with 
various assumptions as to surface energy den- 
sities to be assigned to the shell and compensating 
charges. We begin by assuming that zero surface 
energy needs to be assigned to the shell. 

4. The charges through whose motion the 
energy flow is defined can be regarded as making 
up a net surface charge density o and a net 
surface current density 7; given for a time by 

o=—(1/4r)D-n, (5) 

i,=(c/4r)[HXn], (6) 
where is the unit vector outwardly normal to 
the shell. These equations take care of the transi- 
tion to zero value of normal D, and tangential H. 
Equations (5) and (6) can continue to hold for 
some time without introduction of fresh charge. 
We see this by calculating the flux of current 


out of a closed curve in the surface. This is 
given by 


P 
J -[axn)-as=- - f Cxnyxn-ds 
0 4n Jo 


c 
a - [ ias 
4a 0 
P 
=— - ff cur H-dS 
4n 


= [ fr-as 
Oo | 

=f fr-as- ff —n-as, (7) 
ot 


where J, is the conduction current density. 

The magnetic induction B enters the magnetic 
shell, and is completed as a solenoidal vector 
there by a sheet of induction parallel to the 
surface. It distributes itself there so as to take 
care of the transition of tangential E to zero. 
That is, if d@/dt is the rate of change of the 
surface flux sheet, 


a@/at= —c[EXn]. (8) 


Equation (8) is consistent with B in the volume, 
since the flux of d@/d¢ out of any closed curve in 
the surface is 


OP 

s-| xn|-dS=c f EXn]xn-ds 

0 Ot 0 

—c { B-ds 

—cf f curtz-as 
0B 

ffs. (9) 
ot 


5. According to Eq. (6) the rate of doing work 
by the moving charges is 


c 
[ [zinas—— f [e-Crxnn-as 
4 


-— [ [texi-as. (10) 


We assume that this all comes from energy 
changes and activities within the volume alone. 
Then the Poynting vector P,=(c/4r)[EXH] 
may be regarded as the energy flow density out 
of the region. 

6. The conclusion drawn in the preceding 
paragraph rests on the assumption that no 
changes occur in the energy which may need to 
be regarded as stored in the compensating 
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charges, the shell, or the inner magnetic shell. 
This assumption leads to a self-consistent energy 
flow vector. 

For suppose a large volume is subdivided into 
small volumes, and each volume as well as the 
whole is surrounded by a shell. For each small 
volume we correlate the changes of energy with 
the activity of the adjacent charges in its shell. 
But the activities of the charges adjacent, re- 
spectively, in the same shell to two contiguous 
volumes will be equal and opposite. Hence if the 
activities of the charges in all the shell are 
added up, we will be left only with the activities 
of the charges in the outer shell. Hence the energy 
change in the whole large volume, as defined by 
the activity of the charges in the outer shell, 
is equal to the sum of the energy changes in the 
smaller volumes, as defined by the activities of 
the charges in their respective shells, since no 
energy is assigned to the charges in the shells 
themselves. Thus the energy changes in a given 
total volume as defined by activities of charges in 
shells subdividing it are independent of the par- 
ticular way in which the subdividing shells may 
be constructed. 

7. Instead of moving charges directly along 
lines of current flow as defined by (6), and getting 
work at a rate per cm? given by 


E-i,=(c/4e)[EXH]-n, 


we may move charges in other ways and get 
other expressions for the activity per cm?*. For 
example, if charges are moved about in a small 
circular path in the shell so as to make a current 
of magnitude 7 flowing in a small circle, they will 
produce the same magnetic effect as a magnetic 
double layer of intensity M = 41 per cm? over the 
circle. At the same time, the rate of doing work 
upon the charges will be 


(1/c)(dBn/dt)i= M(dBn/dt) ergs/sec. per cm?, 


of the circle. 

Let M be an arbitrary scalar function of posi- 
tion in space. Let charges be moved about in 
small circles in the shell so as to produce the 
same magnetic effects as a double layer in the 
shell of moment M per cm?. Then H will have a 
discontinuity at the double layer of an amount 
equal to 4rgrad, M. The additional current 
sheet i2 necessary to take care of the transition 
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for tangential H to zero will now be 
c 
tomer wee grad M) Xn]. (11) 
T 


The rate of doing work upon the charges, per 
unit area, is now 


OB, 
E-12+M— 
ot 


c 0B 
= | —[exi]—cLEXerad u}+u—| “Nn, 
An H 


c 
= | “(ExH]~c curl ME}. (12) 
T 


We thus arrive at the alternative Poynting 
vector 


Cc 
P,=—[EX(H—4r grad M)] 


4a 


OB ¢ 
+ M—=—[EXH]-ccurl ME. (13) 
Ot An 


It evidently has the same divergence as P). 

8. We also might move the compensating 
charges about in the shell as follows. Let N be an 
arbitrary scalar function of position. Over the 
inner surface of the shell lay electric double 
layers of moment N per cm*—that is, maintained 
by impressed normal electromotive force of mag- 
nitude N/4xr. Let the moving charges rise up 
through these double layers, with a current 
density J; The rate of this working will be 
— (1/42) NI,-n per unit area. Above these double 
layers the compensating current will be again 
given by Eq. (6). However, tangential E will now 
be changed to E,—1/4m grad, N. The rate of 
working upon the changes per unit area now is 


1 1 
|z-— grad v)-is-—NIen 
4r Ar 


c 1 1 
—_— | (2-5 rad v) xu|-—w1, n 
4n 4a c 
c 1 1 
=“|[ExH]-—ferad nxH}--N1| m 
Ar ar c 


c 1 
=~ | CEx#]-— curl nit} -m, (14) 
T 


T 
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We thus arrive at the alternative Poynting 
vector 


c 1 1 
P;= \| (2- grad w) xu|- vn| 
4r 4n c 


- ~|CExi]- 


4n T 


curl vu}. (15) 


It evidently has the same divergence as P. 

9. It is frequently stated! that the conven- 
tional Poynting vector P; differs from all alter- 
natives such as P.2 and P; in that P; is a function 
of the field vectors alone. This is not true. In 
(13) and (15) we may let M or N be scalar func- 
tions of the field vectors such as a, DE-H, etc., 
where, of course, a and } are dimensional con- 
stants. 

10. We call particular attention to the form 
which P; takes if we set 


—(1/4r)N=¢, (16) 


where ¢ is a scalar potential, and A a correspond- 
ing vector potential of the electromagnetic field. 
That is 


E=-—grad g—10A/cdt, (17) 
B=curl A. (18) 

We then get 
P3=¢1.—(1/44)[(0A/dt)XH]. = (19) 


The first term gJ, is the formula for energy flow 
used extensively by electrical engineers of power 
systems, and where displacement currents are 
only a small part of the total current, may be 
determined directly by wattmeters. We see that 
to make it valid and have the same divergence 
as the Poynting vector, the second term — (1/47) 
X((@A /dt) XH] must be added. 

11. In Section 6 we have seen that by assigning 
all the activity of the compensating charges in 
the shell space to energy changes within the en- 
closed volume and assigning no energy changes 
to the shell and compensating charges them- 
selves, a consistent definition of energy of the 
electromagnetic field may be obtained in that the 
energy assigned to a large region will be the sum 
of the energies assigned to all the sub-regions 


! For emengie, see Jeans, 
fourth edition, p. 519. 


Electricity and Magnetism, 
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into which it may be divided, independently of 
the manner of the sub-division. We have already 
intimated that it is not necessary to be so drastic 
with respect to assignment of energy to the shell 
and compensating charges to get this result. 

We may assign any energy density to the shell 
and compensating charges, provided that upon 
reversal of the direction of the normal to the 
shell the energy density assigned changes sign. 
In this case the argument of Section 6 goes 
through unchanged. At the same time, the 
mathematical expression for the energy density 
assigned to the shell and compensating charges 
may agree with the energy which will be assigned 
to charges and currents within the volume. 

12. For example, in 4 we may assign an energy 
per unit area to the current sheet 7; of Eq. (6) 
of an amount 


1 
T’=-A-i, (20) 
c 
where A is an arbitrary vector function of 


position. If JT is then the appropriate energy 
density for the enclosed volume 


ff fr 

-— f frmas+ f [Stexiny-as 
aS S Agente as 

+f {texas 

-f fo “| xn|+Cex1]| 


We thus arrive at a new alternative Poynting 
vector 


‘dS. (21) 


texi}+—|-4 xn| 
dtLc 


Cc 1 0A 1 0H 
m! | (e+) xt] 4+] x—||. (22) 
4n c ot Cc ot 


Its divergence differs from P,, for 





JOURNAL OF APPLIED PHYSICS 















c2 = fw 


— 


21) 


ing 


22) 
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—div P,= —div P, 


1 @a ar 
-—_(. mm s——A- 1). (23) 
4r ot 


Letting 7, be the energy density determined by 
the conventional Poynting vector Pi, we see 
that 


1 1 
T,=T1T,+-A -I,——H -curl A. (24) 


c 4n 


This equation does assign an energy (1/c)A-4, 
per unit area, to the compensating charges, con- 
sistent with the initial hypothesis of this section. 

13. In particular, we may let A be a vector 
potential of the electromagnetic field defined by 
Eqs. (17) and (18). Then (22) and (23) become 


c | 1 0H 
P,;=—+\ —[grad exH}+-|4 x—|| 
4a c ot 


c 1 0H 
= —— curl gH+ ol x—| (25) 
4n ot 


4r 


So far as its divergence is concerned, Ps is 
equivalent to 


1 0H 
P.=01,4+ 4 x—| (26) 
4 ot 
and (24) becomes 
. 1 
Ts=7T,—-——+-A - I. (27) 
4r ¢ 


14. Again in Section 4, we may assign an 
energy per unit area to the charge sheet o of 
Eq. (5) of an amount 


T” =¢0 (28) 


where ¢ is an arbitrary scalar function of posi- 
tion. If T is the appropriate energy density for 
the enclosed volume 


0 
2 fff 
ot 
0 c 
‘ J [reas f [texan as 
dt 4a 
a ! c 
- - { { -—p-as+ f {exis 
at 4a 4x 
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-ff\-—- Pr aia “texH)| ds, (29) 


We thus arrive at a new alternate Poynting 
vector 
c 1 ¢ 
=—[EXH]—— —(¢D). (30) 
Ar 4n Ot 


Taking the divergence of P;, we find for the 
energy density 


1 
T;=1T,+—D -grad ¢+ ¢p (31) 


Tv 


where p is the charge density. Equation (31) 
assigns to the compensating charges the energy 
go per unit area, consistent with the initial 
hypothesis (28). 

15. Now let ¢ be a scalar potential of the 
electromagnetic field. Then (30) becomes 


Cc 
P;= —— curl gH+¢I, 


T 

1 d¢ 1f0A 

-——p-_|— xx], (32) 
An ot 4nt dt 


where J, is the conduction current density. This 
has the same divergence as 


1 d¢ 1fo0A 
Py=ol.-——D-| xa. . (33) 
4nr at 4rL at 


16. We may illustrate Ps and Ps by a few 
examples. For steady fields 


P.=Ps=o¢l.. (34) 


They, therefore, show energy flowing only where 
there is conduction current, and do not show 
energy flowing in closed paths in empty space as 
P, does. For a stationary magnetic field, pro- 
duced by steady conduction currents in otherwise 
empty space 
H? 1 
T.= ——+-A-I.. (35) 


8x Cc 


The energy assigned to empty space is every- 
where negative. Similarly for a stationary electric 
field produced by stationary charges in other- 
wise empty space 


T3= —E*/8x+ ¢p (36) 
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and again, the energy assigned to empty space is 
negative. 

Consider a plane electromagnetic wave travel- 
ing in the positive z directions and with E and H 
parallel to the x and y axes, respectively, so that 


E.=Eo sin 2nf(t—z Cc), (37) 
H,= Ho sin 2xf(t—2/c), (38) 
Eo= Ho. (39) 


We may choose the scalar potential ¢=0 and for 
the vector potential which will have a component 
along the x axis only 


c 


Eo cos 2nf(1--), 
2af c 


The total current density will all be displacement 
current density and will have a component along 
the x axis only 


A;= (40) 


f 2 
I,=—E, cos 2nf( 1 ). (41) 
2 c 
For the conventional Poynting vector 
c z 
P,=—E,? sin? 2ni( 1 ). (42) 
4a c 
1 Zz 
7 ,=—E,? sin? anf 1 ). (43) 
4a Cc 
For Ps 
c z 
P,.=—E,? cos? anf 1 ). (44) 
4n Cc 
1 2 
T.=—E,? cos? anil _ ). (45) 
4a c 


We see that P, like P; shows the wave as trans- 
porting its energy density with the propagation 


velocity c, but Ps places its energy one quarter 
wave-length behind that of P,. 

Consider now the entirely valid energy flow 
vector 


P,=}(Pi+Ps). (46) 
Then in the example we are considering, 

P 

P,=- —E po’, (47) 
8x 
1 

T3= — Or. (48) 
8 


Thus all periodicity disappears in the energy flow 
and energy density distribution of the electro- 
magnetic wave. However, as soon as the wave 
interacts with matter, the periodicity reappears 
in the energy. For example, suppose the wave is 
absorbed by a current sheet J, at z=0, so that 


c 
[,=——H, sin 2rft. (49) 


4n 
Then Py assigns an energy density to this current 
sheet of an amount (1/2c)A,J, which varies in 


time according to equation 


0/1 Cc c 
“(—A I.) = —F — —F,? sin” 2rft. (50) 
0t\2c 8x 4n 


Comparing with (47), we see then that there will 


be an appearance of energy as heat or otherwise 
given by 


071 
P,——( Au.) 
dt\2c 


Cc 
=—E,? sin? 2xft=P,(s=0). (51) 
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The Use of Secondary Electron Emission to Obtain Trigger or Relay Action 


A. M. SKELLETT 
Bell Telephone Laboratories, Inc., 463 West Street, New York, New York 


(Received June 25, 1942) 


The use of secondary electrons to obtain trigger action similar to that of a thyratron is de- 
scribed. An experimental tube and the necessary circuits by which this action is achieved are dis- 
cussed. This combination gives the features of a triode with a relay or on and off feature, 
resulting in an amplifier, oscillator, modulator, or other vacuum tube device which may be 
turned on or off abruptly at high or low frequencies. In addition, it can be used to replace 
thyratrons in many of their circuits where very low impedance is not necessary and is capable 
of much greater speeds of operation in such applications. 





INTRODUCTION 


URING the past few years techniques have 

been developed which make practical the 
utilization of secondary electrons in commercial 
types of vacuum tubes. These developments have 
been applied primarily to the problem of the 
secondary electron multiplication of current. 
There is described herein another practical use 
of secondary electrons: namely, to obtain elec- 
tronic trigger or relay action. 

The negative resistance characteristic of an 
element that is emitting secondary electrons was 
discovered by A. W. Hull! in 1915 and is fre- 
quently referred to as the dynatron character- 
istic after the name given by him to his tube. 
This was never employed widely, however, be- 
cause of its instability. Recognizing that this 








DEFLECTOR COLLECTOR 
GRID 


1ST ANODE 


Fic. 1. Tube structure. 
1A. W. Hull, Proc. I. R. E. 6, 5 (1918). 
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instability was due to contamination from the 
hot cathode reaching the secondary emitting sur- 
face, and that this contamination travelled in 
straight lines, Jonker and Overbeek? in 1938 de- 
scribed a practical and stable amplifier tube em- 
ploying secondary emission. In their tube the 
secondary emitting surface was shielded from 
the cathode and the primary electrons were con- 
strained to travel in curved paths in order to 
reach it. In the experimental tube described 
herein this same principle is employed to obtain 
stability of operation. 


TUBE STRUCTURE 


The tube structure is shown in sectional view 
in Fig. 1. Surrounding the cylindrical cathode, 
there is a control grid and coaxial cylindrical 
anode of the usual kind but with the modification 
that a hole or slot in the anode allows about 10 
percent of the electrons to pass through. These 
electrons are deflected by the deflector and at- 
tracted by the collector grid so that they follow 
approximately the dotted paths shown. The de- 
flector plate is connected to the cathode and the 
collector grid is maintained at a high positive 
potential. The shape of the deflector plate was 
determined by rubber model studies.* 

The trigger action of the tube depends on the 
properties of the secondary anode which will be 
described further on. This element, bombarded 
by the primary electrons from the cathode, emits 
secondary electrons in a ratio greater than unity 
and these secondaries are drawn off by the col- 


3 J. L. H. Jonker and A. J. W. M. v. Overbeek, W. Eng. 
15, 150 (1938). 

*P. H. J. A. Kleynen, Philips Tech. Rev. [2] 11, 338 
(1937). 
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lector grid which is always positive with respect 
to the secondary anode. 


Properties of the Secondary Anode 


Figure 2 shows a family of current vs. potential 
characteristics of the secondary anode. The differ- 
ent curves correspond to different potentials on 
the collector grid. As the potential is increased 
the current passes through the zero axis twice. 
At each of these two points the number of pri- 
mary electrons is equal to the number of second- 
aries which are drawn off of the surface. If the 
secondary anode is at the higher zero potential 
the external circuit may be entirely disconnected 
and the element will float in a perfectly stable 
manner at this voltage. At the lower zero poten- 
tial, however, the equilibrium is unstable and if 
the external circuit is disconnected the potential 
will “fall” either way. If it falls down the element 
returns to zero and if it falls up it jumps to the 
second zero potential and floats stably there. 

To understand the reason for this stability it 
is only necessary to consider the conditions af- 
fecting the addition or removal of charge from 
the secondary anode. Suppose that initially this 
element is floating at the second zero and that it 
starts to drift lower in potential. The voltage 
difference between it and the collector grid will 
increase and more secondaries will be drawn 


in potential and will thus be brought back to 
equilibrium. Conversely if it starts to drift up in 
voltage the potential difference between it and 
the grid will be decreased and fewer secondaries 
will be drawn across, thereby decreasing its po- 
tential by the accumulation of more negative 
charge. 

Similar reasoning will explain the instability of 
the first zero, except that, whereas for the second 
zero the number of secondaries which are effec- 
tive is that drawn off by the potential difference 
between the collector grid and secondary anode, 
for the first zero all of the secondaries are drawn 
across and the true secondary to primary ratio 
is the important factor. The voltage of the first 
zero is therefore determined primarily by the 
secondary to primary ratio of the surface of the 
secondary anode and is thus approximately the 
same for all values of collector grid potential, but 
the voltage of the second zero is determined 
primarily by the geometry of the collector grid 
and secondary anode structure of the tube. 

The curves of Fig. 2 were taken with the first 
anode current approximately constant at 4 milli- 
amperes. It was adjusted to this value at the 
start of each curve and changed very little. The 
first anode was held at the potential of the col- 
lector grid. 

If the secondary anode is connected to ground 

































































away from it. It will therefore go more positive through a high resistance, electrons will flow 
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Fic. 3. Basic circuit No. 1. 


from the ground through this resistance to it 
because of its positive potential and these elec- 
trons will counteract the effect of an equal num- 
ber of secondaries simply by replacing the charge 
carried away by the latter. Now the effect of 
reducing the secondary to primary ratio is that 
the second zero point or stable floating potential 
is reduced. The secondary anode will now float 
at a potential slightly less than it did when free. 
This new floating potential is determined by the 
intersection of the I-V characteristic with the 
line whose slope is equal to the load resistance. 
For example, the dotted lines R; and R:2 of Fig. 2 
have the slopes of one megohm and two-tenths 
megohm, respectively, and the potentials of the 
points at which they cross the curves are the 
floating potentials for these resistance values. 
Resistances smaller than about 1.4 10° ohms do 
not intersect the curves at all and the critical 
resistance is therefore near this value. It is given 
by the formula 


R.= Em/Im, 
so that in general we have that 
R> En/Inm, (1) 


where E,, and J, are the voltage and current of 
the negative maximum of the characteristic. The 
secondary anode will not float at values of re- 
sistance less than R, for then the number of 
electrons supplied via the resistance will be in 
excess of the number of secondaries needed to 
maintain the effective secondary to primary ratio 
at a value equal to unity and the secondary 
anode will simply charge up negatively and its 
potential will be brought to zero. 


CIRCUIT APPLICATIONS 


These principles may be usefully applied to 
get on and off action in a vacuum tube similar 
to those obtained in a gas tube or thyratron. A 
simple circuit for accomplishing this is shown as 
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Fic. 4. Basic circuit No. 2. 


circuit No. 1 on Fig. 3. The tube is shown in 
schematic form, a break in the first anode being 
indicative of the slot through which the electrons 
pass to the collector grid and finally to the sec- 
ondary anode shown in the upper right quadrant 
of the tube diagram. The values of R; and R: are 
adjusted so that their total adds up to a value 
greater than R., preferably considerably greater. 
A value of from 1 to 5 megohms is recommended. 
Between the cathode and R, is connected a nega- 
tive or C battery of a value sufficient just to cut 
off the flow of electrons from the cathode. The 
grid of the triode structure is connected as shown 
between R; and R». The first anode is connected 
to a positive battery in a conventional manner 
through a load impedance and the collector grid 
is maintained at the highest voltage of the 
battery. 


RI = 1,600,000 OHMS 

R2= 300,000 OHMS 

R3 = 1,850,000 OHMS 
“c" BATTERY = -37 VOLTS 
“B” BATTERY = 150 VOLTS 
NO LOAD IN PLATE CIRCUIT 


eer ad 


PLATE CURRENT IN MILLIAMPERES 
SECONDARY ANODE CURRENT IN MILLIAMPERES 





° 


° 
“2 -10 6 -6 -4 =2 ie) 2 a 6 8 
GRID POTENTIAL IN VOLTS 


Fic. 5. Operating characteristic for circuit No. 2. 
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Fic. 6. Pulse inverter circuit. 


Suppose that the cathode heater is turned on 
and brought to operating temperature. The sec- 
ondary anode and the inner grid will all be at the 
negative C potential and no electrons can flow 
through the tube. If now the potential of the 
secondary anode is momentarily raised to a 
value greater than the first zero potential (Fig. 2) 
the potential of the grid will be raised a few volts 
because of the drop across R; and R: and a small 
electron current will flow through the tube. As 
soon as electrons flow to the secondary anode its 
potential will jump to the second zero voltage 
carrying the grid bias up to the operating point. 
The tube has then triggered on or ‘‘fired’”’ and 
can be utilized as an amplifier through the 
transformers shown. 

It can be triggered off by decreasing momen- 
tarily the potential of the secondary anode to a 
value slightly less than the first zero value or by 
sending a negative pulse into the input or grid 
circuit of sufficient value to cut off momentarily 
the electron flow. The triggering pulses to the 
secondary anode may be applied at A in the 
diagram through a suitable condenser. 

Circuit No. 2 (Fig. 4) is the same as circuit 
No. 1 with the addition of R3. Ri, Re, and R; now 
form a bleeder resistance across the total supply 
battery maintaining the secondary anode at a 
potential slightly in excess of the first zero voltage 
when no electrons are flowing in the tube. R; 
and R: are so adjusted that, as in circuit No. 1, 
the first grid is held just beyond cut-off when the 
tube is in the off condition. This requires a 
slightly greater C battery than for circuit No. 1 
since there will be a positive voltage drop through 
R2 when the tube is off. ' 
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Now if, through the input transformer, a posi- 
tive potential, which may be as little as a frac- 
tion of a volt, is impressed on the grid, enough 
current (less than a microampere) will be passed 
to trigger the tube on; the secondary anode will 
jump to the second zero as in the case of circuit 
No. 1. Circuit No. 2 thus has an advantage over 
circuit No. 1 in that considerably smaller volt- 
ages serve to trigger it on. It may be triggered off 
by applying a negative signal or pulse in the grid 
circuit great enough to cut off momentarily the 
electron stream. 

Figure 5 shows an operating characteristic for 
circuit No. 2. A positive potential of 1.2 volt 
applied between the junction of R; and R: and 
the control grid triggers the tube on and the 
plate current jumps to point Q on the operating 
characteristic, i.e., from zero to 7.75 milliamperes. 
Positive or negative potentials may then be ap- 
plied and the plate current will swing up and 
down as in an ordinary amplifier tube. If a nega- 
tive potential greater than 10.3 volts is applied 
the tube will trigger off as shown. A comparison 
of the shape and slope of this curve with those 
for the triode structure alone (similar to a 
Western Electric 244-A) confirmed the analyti- 
cally determined conclusion that the addition of 
the trigger mechanism has very little effect on 
the operating properties of the triode structure. 
This is mainly because the secondary anode is so 
well shielded from the main electron steam. 

The value of R3 may be obtained by calculating 
the total value of R; plus R:2 as for circuit No. 1 
and then choosing a value for R; such that the 
junction of R; and R; will have a potential with 
no electrons flowing in the tube ten volts or so 
greater than that of the first zero of the charac- 
teristic (Fig. 2). Obviously R; need not be con- 
nected to the highest potential of the power sup- 
ply although in general it probably will be 
simplest to do so for circuit reasons. The off and 
on potentials of the inner grid and the amount of 
its swing may be adjusted by varying the relative 
values of R; and R: and the voltage of the C 
battery. 

The tube in circuit No. 2 is equivalent in most 
respects to a hot cathode thyratron. These differ- 
ences may be noted, however. Since there is no 
deionization time limitation the vacuum tube 
may be made many times faster than the thyra- 
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tron. The vacuum tube may be “‘extinguished”’ 
or triggered off by a negative voltage on the grid 
and while the tube is on, it is not noisy like the 
thyratron and its grid maintains complete control 
over the space current for amplification oscilla- 
tion, detection, modulation, etc. Against these 
advantages should be noted the fact that the 
vacuum tube for the same size cathode has a 
higher internal impedance than the gas tube has, 
or in other words, a relatively larger cathode 
must be used in the vacuum tube to obtain the 
low impedance inherent in the gas tube where 
this characteristic is important. Thus the gas 
tube may be replaced in many of its applications 
by the vacuum tube in circuit No. 2. 

Two typical gas tube circuits, a pulse inverter 
or scale of two counter circuit and a relaxation 
oscillator, are shown as examples in Figs. 6 and 7, 
with the vacuum tube replacing the gas tube. 
The pulse inverter shown makes use of circuit 
No. 1 and the oscillator circuit No. 2. In both of 
these circuits the collector grid is shown con- 
nected to the first anode. It may also be con- 
nected directly to the high side of the battery as 
shown in the circuits above. 

In the pulse inverter, the positive trigger pulses 
fire the two tubes alternately. In the relaxation 
oscillator the condenser C is charged through R, 
until the first anode potential is great enough to 
fire the tube. When the tube fires it discharges 
the condenser and starts a new cycle of operation. 

In the conventional two-tube multivibrator 
circuit the plate of one tube is in phase with the 
grid of the other and is therefore coupled to it 
to give positive feed-back action. In the second- 
ary emission tube shown herein the secondary 
anode is in phase with the control grid so that 
by coupling them together with a condenser a 
one-tube multivibrator may be obtained. This 
circuit is shown on Fig. 8. In operation the grid 
and secondary anode potentials follow wave 
forms similar to those of ‘a multivibrator. The 
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Fic. 7. Relaxation oscillator. 
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Fic. 8. One-tube multivibrator. 


negative grid bias is adjusted to a value that 
would be suitable for amplification. Out-of-phase 
pulses may be obtained across R; or it may be 
shorted out since it is not essential to the multi- 
vibrator action. R2 should be less than R, in order 
that operation may be obtained on the negative 
resistance portion of the characteristic of the 
secondary anode. 


Variations in Tube Design 


Although the tubes described herein have been 
made up with triode structures, it is obvious that 
they may be made up with tetrode or pentode 
structures since the triggering device is external 
to the plate of the amplifier section of the tube. 

For faster operation the opening in the first 
anode or plate may be increased so that half of 
the electron stream is utilized for the triggering 
function. This gives a plate impedance about 
twice as great for the amplifier section. 

If a separate lead is brought out for the deflec- 
tor, this element may be used as an independent 
control of the triggering action. A sufficient 
negative bias on the deflector will prevent elec- 
trons from passing through the slot in the plate. 
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Errata: Electrostatic Electron Lenses with a Minimum of Spherical Aberration 
[J. App. Phys. 13, 49 (1942) ] 


GILBERT N. PLAss 
Research Laboratory of Physics, Harvard University, Cambridge, Massachusetts 


(Received June 10, 1942) 


N the recent paper by the above title, the 
percentages given in the third and fourth 

paragraphs of the ‘“‘note added in proof”’ are in- 
correct. The correct values are as follows: 

For lenses of equal focal lengths and equal 
distances between principal planes, the spherical 
aberration of (14) is 118.5 percent of that of (2); 
and the spherical aberration of (10) is 11.45 
percent of that of (2). This agrees with the 
result that the weak symmetrical lens with 
minimum spherical aberration is (2). The weak 
immersion lens with minimum spherical aberra- 
tion is! 


(A) 


Vi)=A f exp (—32")dz; 


the spherical aberration of this lens is 11.11 
percent of that of (2), only slightly better than 
that of (10). 


For lenses of equal focal lengths and equal 


1 Rebsch and Schneider, Zeits. f. Physik 107, 138 (1937). 


‘‘thicknesses”’ (“‘thickness’’ being defined accord- 
ing to the method discussed in the paper), the 
spherical aberration of (14) is 82.4 percent of 
that of (2); and for (10) it is 40.5 percent of that 
of (2). If the definition is changed, as before, to 
read 1/10 and 9/10 instead of 1/100 and 99/100, 
the spherical aberration of (14) is 94.1 percent 
of that of (2); and for (10) it is 25.0 percent of 
that of (2). 

Using criteria, other than Scherzer’s, for com- 
paring lenses, it is, of course, possible to construct 
“better” lenses than the lens (2). Since it is 
uncertain which criterion would be best, the 
choice between (2) and (14) would have to be 
made experimentally. It was shown in detail in 
the paper that, in principle, these lenses could be 
constructed rather easily. 

I am indebted to Dr. E. G. Ramberg of the 
RCA Manufacturing Company, Inc., and to 
Mr. R. G. Hutter and Dr. L. Marton of Stanford 
University, who have kindly called my attention 
to the errors reported here. 





New Instrument Booklets 








A new stock catalog, No. 18, has recently been issued by 
the Ohmite Manufacturing Company, 4835 Flournoy 
Street, Chicago, Illinois. Rheostats, resistors, tap switches, 
chokes, and attenuators are listed in this 16-page illustrated 
booklet. 


Bulletin No. Z6000, superseding Bulletin No. Z5000, 
provides a condensed listing of the principal measuring 
and control instruments offered by the Wheelco Instru- 
ménts Company, Harrison and Peoria Streets, Chicago, 
Illinois. Notice is given that aluminum casings on all 
instruments have been replaced by steel or cast iron cases, 
for delivery after July 1, 1942. 


Direct current and alternating current indicating instru- 
ments for mounting on switchboards are described in a 
new 12-page catalog, No. 4220, published by the Roller- 
Smith Company, Bethlehem, Pennsylvania. These instru- 
ments include d.c. ammeters and voltmeters and a.c. 
ammeters, voltineters, wattmeters, frequency and power 
factor meters, synchroscopes, rectangular triplex ammeters, 
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and horizontal edgewise triplex ammeters. Round surface 
mounting models and square surface or flush models are 
listed and described. 


Catalog T-625 entitled ‘‘The Homo Method for Tem- 
pering’”’ has just been released by the Leeds & Northrup 
Company, 4934 Stenton Avenue, Philadelphia, Pennsy]l- 
vania. This 28-page, well-illustrated booklet describes 
various uses of the L & N Homo furnace and contains new 
data and new charts. 


An 8-page catalog giving technical data on Emby photo- 
electric cells has been received from the Emby Products 
Company, 1800 West Pico Boulevard, Los Angeles, Cali- 
fornia. A price list is included. 


A new 12-page booklet, giving the history, operating 
characteristics, and details of construction of its mass 
spectrometer, is announced by the Westinghouse Electric 
and Manufacturing Company. Performance data with 
mass spectra for different gases are also given. A copy of 
technical data 85-940 may be secured from Department 
7-N-20, Westinghouse Electric and Manufacturing Com- 
pany, East Pittsburgh, Pennsylvania. 
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